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Abstract 

The effective action of the linear meson model generates the mesonic n-point func- 
tions with all quantum effects included. Based on chiral symmetry and a systematic 
quark mass expansion we derive relations between meson masses and decay constants. 
The model "predicts" values for f v and /„/ which are compatible with observation. 
This involves a large momentum dependent 77—77' mixing angle which is different for 
the on-shell decays of the 77 and the rj . We also present relations for the masses of the 
++ octet. The parameters of the linear meson model are computed and related to 
cubic and quartic couplings among pseudoscalar and scalar mesons. We also discuss 
extensions for vector and axialvector fields. In a good approximation the exchange 
of these fields is responsible for the important nonminimal kinetic terms and the 77-77' 
mixing encountered in the linear meson model. 
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1 Introduction 



In quantum field theory all effects of quantum fluctuations are incorporated in the effec- 
tive action T, the generating functional of one-particle irreducible Green functions. From 
the knowledge of these amplitudes the information about particle masses and decay rates, 
scattering cross sections, etc. can be extracted in a straightforward manner. In particular, 
the effective action for the mesons in the lowest mass pseudoscalar octet contains all in- 
formation on the physics involving only 7r =l= , tt°, K , K°, K and r\. We emphasize that 
all quantum fluctuation effects are already included in the effective action and no further 
integration over fluctuations has to be performed^! Without any further input the effective 
action can be viewed simply as a coherent description of the information gathered by other 
means about scattering amplitudes, decay rates, etc. In a very general context it contains 
already predictive power following from constraints which describe the analyticity properties 
of the momentum dependence of Green functions or general features like convexity. Further- 
more, all exact symmetry relations are automatically embodied in the symmetries of T or 
the related Ward identities. 

Our aim is to find relations among the n-point functions described by T which go beyond 
exact symmetry properties and general constraints. This allows to establish relations among 
physical quantities and to make predictions. (For strong interactions these "predictions" are 
more often "postdictions" , but they permit an understanding of already measured quanti- 
ties.) In the case of meson physics the ultimate goal is a computation of the effective mesonic 
action T from basic QCD, involving as free parameters only a s [Mz) and the current quark 
masses. We will be concerned here with more modest partial answers which follow from a 
few simple assumptions about the general properties of the mesonic effective action. 

A lot of information can be extracted from approximate chiral SUl(3) x SUr(3) sym- 
metry. In the absence of current quark masses for the up, down and strange quark this is 
an exact symmetry of the QCD Lagrangian which is believed to be broken spontaneously 
by the chiral condensate to the vector subgroup SUy(3). Considering the explicit symmetry 
breaking by the quark masses m u , and m s as a small effect and expanding the Green 
functions in powers of these masses gives rise to the very successful chiral perturbation the- 
ory [1], in the context of the nonlinear sigma model for the lowest h octet. In the 
present paper this approach is extended to a linear meson or sigma model including also 
fields for the r/, the lowest lying scalar ++ octet and a scalar singlet. (The latter is often 
called "<r particle", and we use in this work the terms "linear meson model" and "linear 
sigma model" synonymously.) Together with the h octet these fields are combined into a 
complex 3x3 matrix $ which transforms as a linear (3, 3) representation with respect to 
SUl(3) x SUr(3). The corresponding mesons can be interpreted as quark-antiquark bound 
states q L qn [0. In the absence of quark masses spontaneous chiral symmetry breaking arises 
through a nonvanishing vacuum expectation value of the scalar singlet described by the real 
part of Tr$. Nonvanishing quark masses also enforce nonzero expectation values of the 

1 "Effective actions" are aiso often used in a different context where oniy some degrees of freedom are 
integrated out whereas fluctuations of the remaining degrees of freedom still need to be computed. This is, 
e.g., the typical setting of chiral perturbation theory and differs from our approach. 
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diagonal part of the scalar octet. In this context the possible information from approximate 
chiral symmetry breaking is twofold: First, there are a few simple linear relations which 
can be understood easily on the basis of representation theory. A typical example is a Gell- 
Mann-Okubo type mass relation for the particles in the ++ octet. This kind of relation 
could equally well be understood in the context of an extended nonlinear model including the 
rj and the scalar octet. Beyond this, the linear meson model may imply further constraints 
for the free parameters remaining in a nonlinear model to a given order in the quark mass 
expansion. This type of constraint arises typically from nonlinearities in the map from the 
linear to the nonlinear sigma model and is difficult to classify by representation theory. We 
observe that the effective action of the nonlinear sigma model for the pseudoscalar octet is 
completely contained in the effective action for the linear meson model, once restricted to 
the appropriate degrees of freedom. One may therefore hope to extract some information on 
those parameters of the nonlinear sigma model which appear in higher orders in the quark 
mass expansion. 

The predictive power of the linear model is greatly enhanced if approximate chiral sym- 
metry is combined with additional assumptions: 

(i) The derivative expansion assumes for the inverse propagator that the deviation from 

a momentum dependence ~ q 2 + m 2 is only a small effect. This should hold in a range 
of q 2 in the vicinity of the zero at q 2 = —m 2 . It amounts to neglecting terms in the ef- 
fective action which contain more than two derivatives or treating the deviations of the 
inverse propagators from q 2 + m 2 as small corrections in a systematic way. For a deter- 
mination of masses, mixing angles or decay widths only Green functions with on-shell 
external momenta are of interest. Hence, it is natural to expand the proper vertices 
around external momenta corresponding to appropriately chosen "average masses" for 
each SUy(3) multiplet. For many observables this leads to a derivative expansion 
which formally corresponds to an expansion in powers of quark masses. It should be 
stressed, though, that non-analyticities due to multi-particle thresholds clearly restrict 
the range of validity of this expansion. 

(ii) The expansion in the chiral condensate assumes that the typical mass scale rel- 

evant for spontaneous chiral symmetry breaking is small as compared to the typical 
strong interaction mass scales, as, for instance, the string tension or glueball masses. 
The discussion of the relevant scales is somewhat subtle (cf. section ITDp . It is often 



sufficient to assume that the scale of spontaneous chiral symmetry breaking is not large 
as compared to other strong interaction scales. In the present paper we do not exploit 
explicitly this expansion but rather use it in order to establish reasonable ranges for 
some parameters. We emphasize in this context that the present paper makes no poly- 
nomial expansion of the effective action around $ = 0. We will instead expand in the 
difference $ — ($) with ($) the expectation value of $ in the presence of spontaneous 
chiral symmetry breaking and equal quark masses. The latter is justified by the obser- 
vation that a given order in the quark mass expansion only involves a maximal power of 
$ — ($). Within the expansion around ($) the chiral condensate cr ~ Tr ($) appears 
as a free parameter. It should be noted that the validity of a polynomial expansion 
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around $ = would automatically generate a systematic expansion in powers of <To. 
It is, however, not necessary for this purpose. 

Both, the derivative expansion and the expansion in o"o can also be motivated by the 
observation that a classical linear sigma model is in the class of renormalizable theories and 
remains there if it is coupled to quarks. (We neglect here the large-distance nonlocalities 
in the effective quark interactions reflecting confinement induced by the gluonic degrees of 
freedom.) Quantum fluctuations have then the tendency to induce a flow of the effective 
couplings towards the Gaussian fixed point (triviality) in the vicinity of which the deriva- 
tive and polynomial expansions become validQ. Because of nonvanishing meson masses the 
running extends, however, at most over a range somewhere inbetween the GeV scale below 
which the mesons form as quark bound states H and approximately 100 MeV where the pion 
mass acts an an infrared cutoff (with graduation because of the rich scalar mass spectrum). 
Nevertheless, the renormalization effects may be substantial due to the existence of strong 
effective couplings || such that the general form of the effective action may already be in- 
fluenced by the vicinity of the Gaussian fixed point. We do not expect that the derivative 
expansion or the expansion in the chiral condensate <To converge very fast under all circum- 
stances. The associated dimensionless expansion parameters are simply not very small. This 
also holds for the expansion in the strange quark mass in contrast to an expansion in m u and 
rrid- We will discuss these issues in detail and find that the "rate of convergence" depends 
quite significantly on the physical quantity considered. As a general rule, the convergence is 
much better for the flavored mesons than for the non-flavored ones. 

(iii) The leading mixing approximation attributes the dominant deviation from the 
low order results of the quark mass or the derivative expansion to a mixing of states. 
A prominent example is the rj-rj' mixing which indeed turns out to be responsible for 
the comparatively slow convergence of the straightforward quark mass expansion in 
this sector. Another important feature in this context is the "partial Higgs effect" 
which describes the mixing with the h states contained in the divergence d^p'X of 
the axialvector fields. 

The smallest common denominator of all these considerations and the minimal starting 
point for any systematic study of the linear meson model assumes an effective action consist- 
ing of the most general effective potential for <3> and the most general kinetic term involving 
two derivatives. By this we mean that all invariants consistent with SUl{3) x SUr(3) sym- 
metry have to be included which contribute to a given order in the quark mass expansion. 
It is crucial in this respect that the most general kinetic term in the effective Lagrangian is 
not simply Z^Tr d^d^Q. There are other important invariants involving two derivatives 
as, for example, a term 

-X~{ Tr - M $t$) _ <9M$t$) 

+ Tr ($^$ f - d M $$ f ) ($^$ f - } (1.1) 

2 Nonlocal quark interactions related to confinement counteract this tendency for very low momentum 
scales. They influence the mesons only indirectly and are suppressed by a nonvanishing constituent quark 
mass. 
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After chiral symmetry breaking this term induces different wave function renormalizations 

for the pseudoscalar and the scalar octets, a momentum dependent mixing of 77 and rj and 

similar effects which all turn out to be quantitatively important! We emphasize that the 

symmetry breaking effects in the wave function renormalizations Z$ for the various mesons 

(which are described by the kinetic terms) are as important for an understanding of the 

2 

meson mass spectrum as the "unrenormalized mass terms" M i (described by the effective 
potential). With effective inverse propagators ~ Ziq 2 + M i the physical masses are given 

—1/2 

as Mj = MjZj . A study of the mass splitting between the scalar and the pseudoscalar 

— 2 

octet involves the effect of chiral symmetry breaking on M i and Zj. We therefore investigate 
the kinetic terms in the same way as the effective potential. This explains most of the 
differences of our results with earlier investigations [fFfl— 1|14|| where chiral symmetry breaking 
in the kinetic terms was neglected. 

The present paper is devoted to a systematic study of the effective action of the linear 
meson model based on the considerations discussed above. For the pseudoscalar sector we 
take as phenomenological input M w ±, M K ±, M K o, M v r, and f K ±. Using this we compute 
partial decay rates for the 7r°, 77 and rj into two photons as parameterized by f n o, f v and 
f v > as well as other quantities of interest. The perhaps most striking outcome is that M v 
as well as the decay constants f v and are determined essentially as functions of only 
one additional parameter, with a rather weak dependence on the other couplings present in 
the linear sigma model. Fixing this parameter by the measured value M v = 547.5 MeV we 
predict to first order in the quark mass expansion and first order in the derivative expansion 
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1.00 

1.23 

0.91 . (1.2) 



Taking into account the theoretical uncertainties these results are in satisfactory agreement 
with the experimental observations (/W/ w ) ex P = 1.00 ± 0.04, (f v /U) exp = 1.06 ± 0.05 and 
(fr)'/fiz) exp — 0.81 ± 0.02. In view of the lowest order result for vanishing quark masses, 

(f v / fnY°^ = \/3, (fri'/fnY ^ = ^3/8 this is rather remarkable. The values of M v , f v and 
frji for different parameters of the model can be found in sections 0, [| and [13|. One can get 
an idea about the "robustness" of the estimate ( |1.2| ) from the figures and tables of these 
sections. We also discuss the masses of the mesons in the lowest lying ++ octet (section 
|Tl"l) . We find that the scalar partner of the rj has a typical mass of (1300 — 1400) MeV and 
should be associated with the resonance /o(1300) ||19]| . Large mixing effects with two-kaon 
or four-quark states are characteristic for the isotriplet a (980). The resonance a (980) 
may actually be dominantly a two-kaon state and in this case the model suggests a further 
isotriplet resonance with a mass around 1300 MeV. It may be identified with the reported 
resonance ao(1320) fl5|j . 

Four main lines enter our systematic analysis: 
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(1) The relations between the pseudoscalar meson mass differences within a given mul- 
tiplet and the differences in decay constants fx — f w or f K ± — f K o involve the couplings 
between two pseudoscalar octets and one or more scalar octets. Up to a wave function 
renormalization the differences of decay constants correspond to the expectation values (h) 
of the diagonal fields in the scalar octet. Lowest order mass differences follow from the 
cubic couplings ~ Tr(m 2 /i) once the expectation value of h is inserted. (Here m denotes 
the traceless hermitean matrix of pseudoscalar octet fields and h that of the scalar octet.) 
Similarly higher order corrections arise from quartic couplings as Tr(mhmh) and so on. The 
quark mass expansion is closely related to an expansion in powers of the S77y(3)-breaking 
expectation value (h). This mechanism is described in detail in sections ^| and |3|. As a 
byproduct of our analysis one also gains information on the cubic and quartic couplings 
involving pseudoscalars and scalars which are relevant for the decay of a scalar into two 
pseudoscalars etc. 

(2) We formulate the quark mass expansion as a power series in the parameters A u , Ad, 
A s which measure the deviation of the scalar expectation values from their values for zero 
quark masses. In particular, A s — |(A U +Ad) ~ fic — f-n corresponds to the SUy (3)-breakmg 
induced by the mass of the strange quark, and A u — A^ ~ f K ± — f K o measures the amount 
of isospin breaking. In our analysis we actually never need to determine the current quark 
masses. Our "predictions" involve directly the relations between meson masses and decay 
constants. Within the language of a general SUy(3) symmetric model for m, h, etc. the 
essential ingredient is the determination of the cubic and quartic couplings between scalars 
and pseudoscalars as well as the wave function renormalizations from the couplings of the 
linear meson model. This is done in section ^. 

A systematic expression of the pseudoscalar octet mass splitting to order A needs the 
identification of those terms of the effective potential which contribute to this order whereas 
the quark mass corrections to the kinetic term can be neglected. Similarly, an estimate to 
order A 2 (corresponding to second order in the quark mass expansion) involves the effective 
potential contribution to order A 2 , corrections to the kinetic terms to order A and a lowest 
order estimate of the terms involving four derivatives. We find that the apparent convergence 
of the expansion in A is quite satisfactory for the flavored pseudoscalars 7r ± , K ± , K°, K° 
and the 7r°. On the other hand, the formal series in A does not converge very well in the 77-77' 
sector if the singlet mass term generated by the chiral anomaly is, as usual, considered as 
a quantity ~ 0(1). The reason are the relatively large mixing effects which are formally of 
the order A. This is combined with the observation that a zeroth order mass term for the 77' 
(without mixing) is considerably smaller than the physical rf mass and actually not so much 
larger than the zeroth order mass of the 77. We wish to stress that the apparent convergence 
of the A-expansion in the 77-7/ sector improves substantially if one includes systematically 
all effects to a given order in A for all elements of the 2x2 matrix which describes the 
inverse propagator of the 77-77' system. After diagonalization this procedure amounts for the 
mass eigenvalues to a partial resummation of terms which are formally of higher order in A. 
We believe that these convergence properties of the quark mass expansion are quite general: 
The series converges very well only in situations without large effects from mixing of states. 
If mixing is important, a good convergence can only be obtained if the A-corrections are 
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retained for all elements in the relevant matrix. We encounter a very similar situation if we 
want to interpret the a (980) as the isotriplet member of the scalar octet described by h. 

For a computation of differences in decay constants like fx — f n the lowest order in 
the quark mass expansion requires corrections of order A both for the effective potential 
and the kinetic terms. This implies that all relations between meson mass differences and 
decay constants need to lowest nontrivial order also the kinetic terms to order A. The only 
exceptions are the Gell-Mann-Okubo mass relations to linear order in A. They do not involve 
the decay constants and the leading quark mass corrections to the kinetic terms cancel for 
these relations. All other "predictions" of the model involve the quark mass corrections 
to the kinetic term. In order to be able to study separately the quark mass expansion 
and the derivative expansion we have split our systematic exploration of the linear meson 
model to order A into several parts. Sections |] and [| deal with the A-expansion of the 
effective potential. In section |6| we supplement these considerations by a general discussion 
of corrections ~ A and higher derivative contributions to the kinetic terms. The quark mass 
corrections to the kinetic terms are computed within the linear meson model in sect. || and 
the higher derivative contributions are estimated in section [12. 



(3) The expansion in powers of A is a self-contained systematic formalism. Nevertheless, 
large mixing effects sometimes prevent a fast convergence of the series. One of the examples 
encountered in this work is the relatively large momentum dependent off-diagonal element 
in the inverse 77-77' propagator. In turn, the sign and size of this element can be explained 
by the mixing of the pseudoscalar octet m and the singlet p with the states corresponding to 
the divergence of the axialvector fields d^p^. These states have the same quantum numbers 
and the mixing corresponds to the so called "partial Higgs effect". We have performed 
in appendix |B| the corresponding analysis of the vector and axialvector system coupled to 
pseudoscalars and scalars and estimated the contribution from the exchange of d^p^ A to the 
off-diagonal element in the 77-77' inverse propagator. This estimate coincides rather well with 
the value that leads to realistic numbers for M v , f v and f v >\ This is rather encouraging since 
the large O(A) effects find now a natural explanation. Quite generally, our results lead to 
the picture that the A-expansion as well as the derivative expansion converge well once 
a large enough basis of states is included. Integrating out such states, however, can lead 
to large coefficients in the formal A-expansion for the remaining states and therefore to a 
slow convergence. Keeping the additional states or integrating them out without a further 
truncation of the series to given order in A is equivalent to a resummation of higher order 
terms from the point of view of the formal A-expansion. This improves the convergence 
substantially. We therefore have given a general discussion of mixing effects in appendix 
where we also identify the most prominent mixings relevant in our context. 

(4) The effective action T includes all quantum fluctuations. If T is known no further loop 
calculations are necessary. Nevertheless, it is often useful to estimate the contributions of 
fluctuations for certain invariants contained in T. An example are the contributions to higher 
derivative terms which reflect the deviation of an inverse propagator G~ 1 (q 2 ) from the leading 
momentum dependence Zq 2 + M . Recently, a reformulation of perturbation theory was 
based on effective vertices and propagators instead of the classical ones |16| . Here, this means 
that the contributions of quantum fluctuations to the difference G~ l (q 2 ) — (Zq 2 + M ) can be 
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computed in terms of Z, M 2 and appropriate effective cubic couplings 7*. Since the effective 
vertices have been determined by our analysis we can give in sect. a quantitative estimate 
of the contributions from scalar and pseudoscalar fluctuations to the higher derivative terms. 
They turn out to be substantially smaller than those arising from mixing effects with other 
states. 

Besides its main purpose of a systematic discussion of meson masses and decay constants 
within the framework of the linear meson model this work also constitutes the basis for 
several interesting developments: Recently, the parameters of the two-flavor meson model 
were estimated from a solution of nonperturbative flow equations ||. The resulting values for 
f w and the chiral condensate (ipipj encourage an extension of this approach to three flavors. 
In this case the nonvanishing mass of the strange quark plays an important quantitative 
role and needs to be included. The results of such a future estimate of the parameters of 
the linear meson model can then be compared with the values infered in the present work 
from a phenomenological analysis. The present results also give an idea which invariants are 
important and should be retained in the necessary truncation of the exact nonperturbative 
flow equation. 

A determination of the parameters of the linear sigma model in the zero temperature 
ground state is the starting point for any analysis of the high temperature behavior of this 
model. This is relevant for the QCD phase transition in the early universe and for ongoing 
and future heavy ion collision experiments. Within the linear meson model the temperature 
dependence can be studied along similar lines as in |TTJ. On a crude qualitative level one 



may also use a mean field approximation fTj-jnj. The temperature dependence of meson 
masses - in particular the vector mesons - may be experimentally observable in heavy ion 
collisions. We should point out, however, that such a study remains reliable only as long 
as 2ttT is smaller than the compositeness scale above which a mesonic description of strong 
interaction physics becomes invalid. 

Finally, we have not exploited in this work the information we have gained concerning 
cubic and quartic meson couplings. Especially for the extension of our model to include 
vector and axialvector mesons (appendix ||) these couplings determine the decay rates and 
branching ratios. They therefore contain a whole lot of additional "predictions" that can be 
compared with experiment. We plan to come back to all these issues in future work. 



2 Meson masses to linear order 

In the limit of vanishing quark masses the chiral symmetry SUl(3) x SUr(3) is spontaneously 
broken to the vector-like SUy{3) subgroup. With respect to this symmetry breaking the 
18 real fields contained in the complex (3, 3) representation decompose into a pseudoscalar 
octet plus singlet and a scalar octet plus singlet. The pseudoscalar octet is associated with 
the eight (pseudo-)Goldstone bosons 7r =t , n°, K ± , K°, K , 77 and the pseudoscalar singlet 
with the rf particle. The scalars comprise the isotriplet a , the strange scalars K^, Kq°, 
K as well as two states with /o quantum numbers, one of them being the "cr-particle" . 
Chiral symmetry breaking is induced by a nonvanishing expectation value 00 of the scalar 
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singlet. In presence of current quark masses the S77y(3) symmetry is explicitly broken. This 
is reflected by nonzero expectation values of the flavor neutral fields within the scalar octet. 

This and the next section are devoted to a general analysis of the effects of S77y(3)- 
breaking without explicit reference to the linear sigma model and spontaneously broken 
chiral symmetry. These concepts will only be introduced in sect. |j. In this way we can 
understand the new relations emerging from the linear sigma model in the context of a more 
general setting involving only S77y(3) symmetry. 

The octet of pseudoscalar mesons is described by a hermitian traceless 3x3 matrix m, 
whereas the scalar octet is similarly denoted by h. In addition, we have the pseudoscalar 
and scalar SU(3) singlets p and s. These fields are normalized with standard kinetic terms 

£ k in = ~ Tr d^md^m + ^Tr d» hd^h + ^pd^p + sd^s . (2.1) 

With SU(3) generators X z obeying Tr X y \ z = 2S yz we may write m = m z X z , h = h z \ z with 
m z , h z real. For vanishing quark masses the pseudoscalar octet corresponds to the massless 
Goldstone bosons. The most general mass term consistent with 577(3) symmetry, charge 
conjugation C and parityf] V reads (we will work in Euclidean space time throughout this 
paper) 

£ mo = -m 2 m Trm 2 + -m 2 h Tr h 2 + -m 2 p p 2 + -m 2 s (s - s ) 2 . (2.2) 

Small nonvanishing quark masses can be described by a linear perturbation or external source 
for the scalar mesons ^ 

Cj = --Tr j h h - j s s (2.3) 

with jh diagonal and traceless. Both jh and j s are linear in the quark masses with[] 

/ 

jh + C s j s ~ 



m d I . (2.4) 



V m s 

This leads to a shift in the expectation value of the scalar fields 

(s) = u 

(h) = w\ 3 — v^3fA 8 



;2.s) 



There are direct relations between the parameters v and w and the differences of meson 
decay constants which are explained in appendix ^ 



^r) 2 v = l(2A s -A u -A d ) = ±(f K± + f K0 -2f„ 



i 

Z m \ 2 



z. 



) 2 w = (A„ - A d ) = f K ± — f 



(2.6) 

K° ■ 



3 Charge conjugation corresponds to a transposition of both matrices m and h whereas parity amounts to 
to — > —to, h — > h, p — > —p, s — > s. 



4 The normalization c s can be inferred from sect. ^. 
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Here A u is proportional to the difference of the (uu) condensate for nonvanishing and van- 
ishing quark mass and similarly for A^ and A s . The quantities are related to the meson 
decay constants f\ by SU (3) breaking wave function renormalizations which will be discussed 
in section Until then we take into account the effects of nonvanishing quark masses only in 
the lowest order in a derivative expansion. In this approximation the kinetic term is given by 
Q2.1| ) and we may identify f i with /j. The wave function renormalization constants Z m and 
Zh are induced by different kinetic terms for the pseudoscalar and scalar octets within the 
linear sigma model and will be explained in sect. [|. The isospin violating expectation value w 
is very small and we can use the experimental values f n = f w ± = 92.4 MeV, /#± = 113 MeV 
to estimate 

Z \h 

"' ' v = 13.7 MeV . (2.7) 



Z 



We will see later, ( 6.21] ), that the relation ( |2.7| ) is modified by quark mass corrections to the 



kinetic terms leading to f i ^ fi. As a result the value (|2.7|) will be shifted to 23.3 MeV. 

At this stage u is not yet fixed, since we still have to specify s or, equivalently, the 
meaning of s = 0. A possible choice would be that s = denotes the minimum of the 
potential in the absence of quark masses, i.e. so = 0. For this choice vanishes and 



(s) = u = y 2/3 (Au + A d + A s ) (Z s /Z m )^ . It will, however, prove to be more convenient to 
choose the origin of s corresponding to the minimum of the singlet field in the presence of 
quark masses. In this case one has s = — j s / m s an d 

(s) = u = , m 2 m > (2.8) 

such that ~ O(A). We note that more generally the choice of s fi xes the point around 
which the potential U(m, h,p, s) is expanded and therefore rr? m as well as the values of all 
other parameters in a polynomial expansion of U depend on s - F° r the choice u = we 
will use m^j, v and w instead of the current quark masses m u , m d and m s as the parameters 
characterizing the explicit chiral symmetry breaking. Their relation to the quark masses 
involves the solution of the field equations for s and h in presence of the sources ( |2.3|) . With 
the choice u = our expansion parameter is given by (Z m / Z h ) l l 2 v ~ \{f k ~ fn)- The quark 
mass expansion turns into a Taylor expansion in the small parameter (f K — /„.)/(/# + 
/„.) ~ 0.1(0.2). Here the numbers in brackets include the effects of quark mass corrections 
to the kinetic terms. Correspondingly, the small parameter for isospin violating effects is 
(f K ± — f k°) / (f k± + f k°) — ~ 2 • 10 _3 (— 3 • 10~ 3 ) (see below). This implies that linear isospin 
breaking effects give corrections of the same order of magnitude as quadratic quark mass 
effects. As a rough estimate one expects that the squared meson masses can be computed 
to linear order in the quark masses with an accuracy of 10-20 percent. This corresponds to 
a typical uncertainty for the masses of the pseudoscalar octet of around (30-60) MeV. To 
quadratic order in A this error is expected to decrease to a few MeV. 

We want to investigate the dependence of the pseudoscalar meson masses on the param- 
eters m^, v and w. We will work in this and the next three sections with a kinetic term of 



the form (|2.1|) and discuss effects from modifications of the kinetic term in sect. |6|. To linear 
order in the quark masses we need the contributions linear in mL, v and w. This involves 
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the cubic couplings of two pseudoscalars and h or s. Their most general form consistent 
with SU(3) symmetry, charge conjugation and parity reads 

1111 
£ 3 = ^7i-s Tr m 2 + -72 Tr m 2 h + -73P ^ r m ^ l^ iS P • (2-9) 

We note that because of (s) = only and the couplings 72, 73 contribute to the pseu- 
doscalar mass matrix. We find for the masses of the off-diagonal or flavored mesons tt , 
and K°, K° 

2 X 

m k± = rn 2 m + -j 2 (v + w) (2.10) 



2 1 

M K0 = m 2 m + -72 (v-w) . 

Again, the bars indicate that the quantities Mi differ from the physical meson masses Mi by 
SU(3) breaking wave function renormalization effects to be discussed in section |6]. As long 
as we restrict the discussion to the kinetic term Q2.1|) (sections we can identify Mi with 
Mi. Eq. (|2.10j) allows one to express couplings in terms of meson masses 



m 2 



i (mI ± + m 2 k± + m 



•2 

K° 



l2W = M 2 K± -M 2 K o (2.11) 



l2 v = -(M 2 K± + M 2 K0 -2M 2 n± 



and to estimate the isospin violation to leading order in the quark masses 

z m y fz m y m 2 k0 -m 2 k ± 



f K0 -f K± = -(^ w = 3 ^ v ^0.47 MeV (2.12) 

which is in good agreement with the result fj^o — fx± — 0.45 MeV from chiral perturbation 
theory 0. Here we have used electromagnetically corrected masses M^± = 135.1 MeV, 
M K ± = 492.4 MeV, M K o = 497.7 MeV. Including quark mass corrections to the kinetic 
terms these values are shifted to (Z m /Z h ) l l 2 w ~ -0.67MeV, (07231) , and f K o ~ 113.28 MeV, 
( |6.24| ). We will occasionally also use the isospin means f K = |(/ K ± + f K o) and M 2 K = 
l(M 2 K± +M 2 K0 ). 

To obtain the masses of the neutral pseudoscalars n°, r\ and 7/ we have to diagonalize 
the general mass term 

- Km 2 + M 2 8 m 2 8 + M 2 p 2 ) + M 2 38 m 3 m 8 + M 2 p m 3 p + M 8 2 p m 8 p (2.13) 
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which also involves m 2 and 73 according to 



Ml 



ni„ 



M 2 3p 
Ml 



Mi 



wL + 7 2 w 



- (2M 2 K± + 2M 2 K0 - M 2 n± 



(2.14) 



m 2 m - l2 v 



m; 8 = 



V3 



72 w 



= m\ 

7T- 1 - 

-L(m 2 k± -m 2 



K° 



We recover the usual lowest order relations of chiral perturbation theory or the "eightfold 

2 2 2 

way" for M 8 , M 3 and M 38 [Tj| Q]. We will neglect the isospin violating mixings of 777,3 which 
only contribute to order w 2 to the mass eigenvalues. Diagonalization of the (m§,p) sector 
yields the mass eigenstates 

77 = 777 8 cos B p — p sin P 
7/ = p cos Op + 7Ti 8 sin Op . 
The masses of 77 and 77' and the octet-singlet mixing angle P are given by 



(2.15) 



M , + M 



M 2 + Ml 



M 2 



Mi 



2 \ 2 

M 



+ AMt p 



Ml 



tan 6, 



M 2 + 



(Ml 



M 2 
p 



(2.16) 



Contrary to the meson masses the mixing angle p will receive additional contributions 
~ O(A) from modifications of the kinetic terms discussed in section ^ 

It will later be our aim to find relations among the parameters ?77p, 777^, 72, 73. For the 
moment we only notice that the couplings 72 and 73 are directly related to the partial decay 
width of the scalar octet into two mesons belonging to the pseudoscalar octet or singlet. In 
particular, 72 can be extracted from 



72 



and ( 2.11 ) and one finds 
11040 MeV . 



(2.17) 



a/2 



72 



4942 MeV once quark mass corrections to the kinetic 

J)- To 

nent O by 

^75sTr/i 2 + ^7 6 Tr/i 3 + ^ 77S 3 (2.18) 
4 6 



This value changes to (Z^/Z, 
terms are included. 

We next turn to the scalar masses. To zeroth order they are given by ml in 
linear order in A we have to supplement (|2.9|) by 

A£ 3 
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leading to 

M 2 a ± = m 2 h -2 l6 v 

M 2 K ± = ml + l6 {v + w) (2.19) 

M K * = m\ + 7e (f - w) . 

In the following we will neglect isospin violation in the scalar sector and use the isospin 
means Mk* and M ao . The diagonal part of the mass matrix for the flavor neutral scalars 
reads 

M 2 S = m 2 s 

M 2 a3 = m 2 h -2 l6 v (2.20) 

M % = m l + 2 lev 
whereas the off-diagonal elements are given by 

Ml = l5 w 



M l = -v 7 ^ ( 221 ) 

We note that 77 does not enter these O(A) expressions for the pseudoscalar meson masses. 
In particular, we find the interesting relation 



M% = - (AM 2 KS - M 2 a ) (2.22) 
which is the Gell-Mann-Okubo formula in the scalar sector. 



3 Pseudoscalar meson masses to quadratic order 

Before proceeding to an estimate of the various couplings we will analyze in this section the 
general structure of the pseudoscalar meson masses to quadratic order in v and w, assuming 
a kinetic term of the form ( |2.1| ). We stress, however, that the modifications of the kinetic 
terms for nonvanishing quark masses are of the same order as the effects discussed in this 
section. They are postponed until sect. || only for the sake of a simple presentation and for 
separating clearly different orders in the derivative expansion. Mass corrections quadratic in 
v and w involve the quartic couplings for two pseudoscalars and two scalars. Their general 
form is given by 

£4 = -5is 2 Ty m 2 + -62s Tr m 2 h + -dsspTr mh + -54S 2 p 2 + -55p 2 Tr h 2 

+ -5qpTt mh 2 H — 67 Tr h 2 Tr m 2 H — 5s Tr m 2 h 2 -\ — <5g Tr (mh) 2 -\ — <5io (Tr mh) 2 . 

11 S (3.1) 
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We note that the couplings 81, . . . , 84 do not contribute to the pseudoscalar masses. The term 
~ 5 5 adds effectively to m 2 an additional piece S 5 (3v 2 + w 2 ) whereas m 2 ^ is supplemented by 
S 7 (3v 2 + w 2 ). The quartic terms contribute to the flavored meson masses 

AMl ± = 5 7 (3v 2 + w 2 ) + 28 8 (v 2 + w 2 ) + 25 9 (v 2 - w 2 ) 

AM 2 K± = 5 7 {3v 2 + w 2 ) + 5 8 {5v 2 - 2vw + w 2 ) - 45 9 {v 2 - vw) (3.2) 

AM 2 K o = 5 7 (3v 2 + w 2 ) + 5 s (5v 2 + 2vw + w 2 ) - 48 9 (v 2 + vw) 

whereas the mass matrix of the neutral pseudoscalars receives corrections 
AM 2 = 5 7 {3v 2 + w 2 ) + 2(5 s + 5 9 ){v 2 + w 2 ) + 5 10 w 2 
AMg = 5 7 (3v 2 + w 2 ) + l(5 8 + 5 9 )(9v 2 + w 2 ) + 35 w v 2 
AM 2 = 5 5 (3v 2 + w 2 ) 

—2 1 (3.3) 

AM 38 = - — (45 8 + 45 9 + 35 w )vw 

AM 2 p = -25 6 vw 

AM 2 8p = -±=5,{3v 2 - w 2 ) . 

The lowest order relations between the off-diagonal and diagonal mesons are only modified 
by the 5 9 and 5±o terms 

M\ = ~Ml± + A5 9 w 2 + 5 w w 2 



Mg = i (2M 2 K± + 2M 2 K0 - M'±) + 125 9 v 2 + 35 w v 2 . ^'^ 



The combination 5 9 + Siq/4 is therefore particularly interesting for the mass relations. 

If we neglect isospin violating contributions of order w 2 to the mass eigenvalues and treat 
the flavored meson masses M n ±, Mr-±, M^-o as three input parameters we can, in principle, 
predict M v and M v i as well as the octet-singlet mixing angle 9 p . This requires information 
about the coupling 5 9 , the off-diagonal mass term 

MI p = -V3v(k + 6 6 v) (3.5) 

as well as the singlet mass term 

M 2 = m 2 p + 35 5 v 2 . (3.6) 

It will be the aim of the next section to discuss these quantities in the context of a linear 
sigma model. In addition, we can use the relations ( [2.10| ) and ( |3.2| ) to obtain information 
about v, w and therefore about the differences of decay constants ([2.6|). Without isospin 
violating effects we have at this stage four unknown observables (M v , M v >, 9 P , fx± — fir) for 
which we want to find relations. 
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Finally we note that the number of parameters can be reduced by absorbing 85 and 87 into 
the definitions of m 2 and m^, respectively. This will be done in section [5] in the framework 
of the linear cr-model. For w = also 8q can be absorbed into 73. The four observables 
depend on the "couplings" + 38 5 v 2 , 73 + 8 e v, 8 S and 8$ + 81O/ '4. 



4 Linear meson model 

In this section we will compute the couplings of sections in the context of the linear 
meson model (often also called linear sigma model). The fields m, p, h and s are all contained 
in a complex 3x3 matrix $ which transforms as (3, 3) with respect to the chiral flavor group 
SU L (3) x SU R (3) 

<$>^U R m{- U R eSU R (3), U L eSU L {3) (4.1) 

and carries nonvanishing axial charge. Including up to two derivatives the effective action of 
the linear cx-model can be written as a sum of a potential and a kinetic term plus a source 
term 

r[$] = J d 4 x(C/ + £ kin + £ i ) . (4.2) 

As a consequence of the invariance under SUl(3) x SU r {3) symmetry and the discrete 
transformations^ V ($ — > *$>*) and C ($ — > $ T ) the potential is a function of the four 
independent invariants || 

p = Tr$ f $ 

r 2 = ^Tr($t$_^ 2 = 3 Tr ( $ t $ ) 2 _l p 2 



r 3 = Tr(V$-ip) =T r($t$) 3 -^r 2 p-ip 3 
£ = det $ + det . 



(4.3) 



With respect to the vector-like SUv{3) symmetry we may decompose 

* = ?0 + -L(» p + ^ + * i + -L,,) (4.4) 
with traceless hermitian 3x3 matrices $ p , $ s and real singlets (<t is a real positive constant) 

1 = [Tr ($ + $ f ) - 6a 

(4.5) 

TV f cf> _ dit^ 

v/6 



X P = -^=Tr($-$t 



5 More precisely, the transformation <f> — * $^ corresponds to left-right symmetry which is closely related 
to the parity reflection V. 
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The kinetic term involving two derivatives consistent with SUl{3) x SUr(3) symmetry, 
C and V readsf] 



-X~{ Tr ($ f ^$ - ^$ f $) - ^$ f $) 

+ Tr ($d M $ f - ^$$ f ) ($^$ f - ^$$ f ) } 
^X+{ Tr ($ f ^$ + + ^$ f $) (4.6) 

+ Tr ($9 M $ f + ($0"$ f + 0"$$*) } 

Tr | (fi^$t$^$t$ + $t^$$t^$) f $t$ _ I Tr 



Here Z^,, V^, etc. are functions of the four independent scalar SUl(3) x SUr(3) invariants 
( |4.3|) and the dots stand for other independent terms which are not relevant for our purposes, 
as for example d^T^d^- (See sect. |8] for the precise meaning of this statement.) We note 
that the additional pseudoscalar invariant 



CO 



(det$-det$ f ) (4.7) 



will always appear with even powers, since there is no other parity-odd invariant. Since its 
square can be expressed in terms of 



uj 2 + i 2 = 4det($ f $) = 4exp{Trln($ t $)} . (4.S 



it can only appear as an independent quantity in combination with derivatives as in 
Evaluating (|4.6|) for a configuration $ = diag(cfo) the structure of the kinetic term for 
fluctuations of $ around $o simplifies to the form 

£ kin = l -Z m Tr d»%d,% + l -Z h Tr d^ s d^ s + l -Z p d» Xp d»X P + \z s d» Xs d»Xs (4.9) 

where the normalization is adapted such that to lowest order in Wo one has Z m = Z% = Z s = 
Z p = Ztp. The fields $ p , Xpi Xs have the same transformation properties with respect to 
SU (3) and parity as m, p, h, s, respectively. 

So far cfo in (4.4) has not been specified. If we want to use an expansion around (s) = 
we should identify it with the value of | (Tt $ + Tr $^ at the potential minimum in presence 
of the quark masses, i.e. 

°o = \ K + a d + a s ) Z^ = 1(7,+ Jk± + Jk°) ■ (4-10) 



6 By partial integration we bring all contributions to the kinetic term into the form where the two deriva- 
tives act on different fields. The invariance of the last term follows from V a a'Vbb'V CC 'e a 'b'c' = tabc for arbitrary 
V e 517(3). 



15 



This corresponds to our choice u = in section |] and leads to the identification (cf. Q2.1J) ) 



,, - (2Z m ) *m , $ s = (2Z h ) *h 
X P = Z p 2 p , Xs 



Z s ~ 2 s . 



(4.11) 



Corrections to the kinetic terms which are linear in the quark masses involve cubic terms 
like Tt ^gd^^pd^p etc. and will be discussed in section RL 

For the configuration $ = diag(<T ) the invariants ( |4.3| ) take on the values 



Po 
Co 
r-2 



3a 2 

r 3 = 0. 



(4.12) 



We next decompose the invariants p — po, r 2 , r 3 and £ — £ into the irreducible representations 
$ p , $ s , Xp an d Xs- Using occasionally the shorthand notation 



A s 



<?o + Xs 



(4.13) 



one finds 
P~ Po 

r 2 



I'M 2 
2 p 



woXs 



(4.14) 



+ 



I Tr $ 2 $ 2 - 1 Tr ($ s $ p ) 2 - 1 Tr $ 2 Tr $ 2 + x s Xp Tr $ s $ p 



4 



4 



(4.15) 



73 



v/3 



>/3 



>/3 



+ ^-Xs Tr <&f + Tr ^ + ^- Xp Tr <^ + Tr 



Tr$ 2 + Tr$ 2 r 2 



72 



Tr $: + Tr $; 



+ 
+ 



-Tr 

8 

1 

4^3 



V3 



(Xs^s ~l~ Xp^p) 



X s Tr ($ 2 $ p $ s $ p - + Xp Tr ($ 2 $ s $ p $ s - & p $\ 



(4.16) 



-I ^Tr $f$ 2 + $ 2 $ p $ 2 $ p - 2^^$,$, 



Including terms quartic in m, h, s, p the invariant r 3 only contributes 

r 3 (4) = 2^2^ Tr $ 3 S + a 2 [3 Tr + $ 2 $J) - Tr $ 2 (Tr $ 2 + Tr $ 

+2v / 3(x s Tr$3 + Xp Tr$ 2 $ p 



(4.17) 
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Using the fact that for an arbitrary traceless 3x3 matrix T 



one also obtains 



det(a + T) = ^ Tr T 3 - Tr T 2 + a 3 (4.18) 



t-t* = ^Tr$ 3 -i=Tr<|. s $ 2 -^o(Tr<|. 2 -Tr$j) 



- ^Xs(Tr®s-Tr®l) + ^X P Tr$s%-voX 2 p -^XsX 2 p (4.19) 

+ V&alxs + ooxl + ■ 

For the choice (|4.10| ) one has u = (c.f. section |2|) and we denote v and w collectively by A. 
One observes that the expectation values of the invariants in the presence of quark masses 
obey 

(p-po)-O(A 2 ) , <£-£o>~0(A 2 ) 

( 4 -2°) 

(r 2 )~0(A 2 ) , (r 3 )~0(A 3 ). 

The mass squared matrix for the pseudoscalar mesons is obtained from the potential of the 
linear cr-model by taking the second derivatives with respect to the parity-odd representa- 
tions $ p and Xpt which we will collectively denote by <p~. Since the potential is parity-even, 
all VEVs of single derivatives of the invariants ([4.3|) vanish. We note that 




x dip~dip~ 

' <9 2 T 3 

dip- dip- 



0(1) 

(4.21) 

0{A 2 ) . 



The most general SUl(3) x SUr(3) symmetric potential can be expanded in a Taylor 
series around p = p Q , £ = £ > r 2 = and r 3 = 0. We now see that a determination of the 
pseudoscalar masses up to quadratic order in the quark masses requires 

U = m 2 g (p- p ) - ~F [£ - £ - ^o(p - Po)] 

+ ^Ai (p - po) 2 + ^A 2 r 2 + ^A 3 r 3 (4.22) 

+ ^i(p-Po)(e-eo) + ^ 2 (p-po)r 2 + ^ 3 (e-eo)r 2 + ^ 4 (e-eo) 2 + .... 

Here the potential has been normalized such that it vanishes for the configuration $ = 
diag(ao). We see that to zeroth order in the quark masses the only contributions to the 
pseudoscalar mass matrix arise from V and m 2 . To linear order we obtain corrections from 
A 2 and V whereas to quadratic order the other Aj and the fi { enter. We add to the effective 
action a source term 

£ J = -iTr(Vj+/<l>) (4.23) 
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which is linear in the real quark mass matrix M q = diag(m u , m^, m s ) 

3 = i f = a-qMq . 
We denote the singlet part of the source by 



and require 



_d_ 



3 > I *=CT 



Our choice ( [4.10 ) for o"o therefore implies 



V6(i m 2 = j s Zi 



(4.24) 



(4.25) 



(4.26) 



(4.27) 



and the mass term m 2 g is linear in the quark masses. 

Comparing ( 4.22 ) with ( |2.2| ), ( |2.9| ) and ( |3.ip we can now determine the various couplings 
of the last two sections in terms of those of the linear sigma model. For the pseudoscalar 
mass terms of (12.21) we find 



m 



2 2 ry-\ 

— m g Zj m 



m2 g + o ua °) Z p 1 



(4.28) 



The cubic couplings of 

7i 

72 
73 
74 



contributing to 0(A) read 
- (3a A 2 + 77J Z h 2 Z m l 



_i _i _i 

— — I C0A2 — — V ) Z v 2 Z h 2 Zrn 2 



(_ -r 1_ 



(4.29) 
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whereas for the quartic couplings of (|3.1|) we obtain 
1^ 5 , 



S3 = (^2 + ^0^ + 3^ + 3^3+^)^^;^-!^ 

Sa = {^i-a^-Aa^Z^Z; 1 

Sb = QAi-^+^ + ^-S^)^- 1 ^- 1 

s 6 = ^(\2 + ^l\ 3 )zpz h 1 z m ^ 



(4.30) 



m 



<*8 = o ( A 2 + 2a A 3) Z h Z, 



111 



S9 — -t^^Z^Z,^ 



3 

S w = 0. 



We next turn to the scalars for which we wish to relate the couplings mi?, m\, 75, 76 to 
those of the linear sigma model. For this purpose we note that, contrary to the case of the 
pseudoscalar mesons, the parity-even scalar meson fields, collectively denoted by (p + , may 
also appear to odd powers in the invariants p, £, r 2 , r 3 . We therefore need in addition to 

l 8 »\ 0(1) , (^)~0(A) 



d(D + / \ <%9 H , 

(4-31) 

^ \ /r\f -\ \ I dT 3 \ ,n/ A 2\ 



and 



\ <9y2+ / \ dip 

9V ^ <XU . (^)~o ( v, 



\ d(p + dip + / \ d<p + d(p^ , 

\d<p+d<p+/ y ' ' \dip+dip+ 1 1 ; 

Hence, the expansion ( |4.22|) of the potential contains exactly those terms required to obtain 
the scalar masses to linear order in A. Furthermore, we see that to zeroth order in A only 
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nig, v, Ai, A2, Pi and /3 4 contribute. We find in particular 



and 



75 
76 



i 2 s = m 2 g Z; 1 + 6w (w X 1 + alf3 1 + a 3 f3 4 -^v)z; 1 



m ft = m Z h + cr (3o- A 2 + H Z 



4 



77 = V6 (|^oAi + |0"o^! + 3<rg/3 4 - — 1/ ) Z 



2 

s 



(4.33) 



t (9ao A 2 - V + 12^ A 3 ) (4.34) 



5 Parameters of the linear meson model 



In this section we will give a first estimate of the values of the parameters of the linear 
a-model. It is based on an expansion in powers of A to lowest order in the derivative 
expansion[]. Comparison of the estimates of this section with those of the following ones will 
allow us to evaluate the quantitative influence of quark mass corrections to the kinetic terms 
(which are neglected here). The results of this section can also be compared with earlier 
work 0-|l4| by setting Z^ = Z p = Z m . Later we will see, however, that Zh/Z m deviates 



substantially from one. We observe that the couplings Ai, (3 1: /3 2 , (3 3 and /3 4 influence the 
meson masses only through m 2 , £5 and 67 whereas A3 appears in addition in 5 6 and <5g. 
The couplings £5 and £7 modify the relation between the neutral and flavored pseudoscalar 
masses, (|3.3| ) and ( |3.4|) , through the term 365V 2 in M 2 and the term 367V 2 in M 3 , M 8 and 
the flavored pseudoscalar masses. (We neglect here corrections ~ w 2 .) A redefinition of 
couplings 

" + (35 7 + A5 8 - 25 9 



m 2 g 



v Z m 



A 



V + 



A, 



25$Z P — 



25 7 



25 5 Z P + 



3 8 

07 + 7^8 



3 J 

-h 

3 - 



(5.1) 



a 2 



7 The systematic ordering of the derivative expansion is ambiguous to lowest order since a minimal kinetic 
term must always be included. For our purpose we consider to lowest order the kinetic term (4.9). The first 
order comprises the most general terms with up to two derivatives, the second order includes four derivatives 
and so on. 
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absorbs this correction in the lowest order masses m m , m 2 , m\, whereas the corresponding 
shifts in the 7's only contribute to cubic order in A. In terms of these couplings one has 

mI ± = ^ 2 z~ l _i (3^0 + 1/) zj (J K±+ J K0 -2J W ) 



6 

- (a' 2 + A 3 ao) Z m 2 (J K± + f R0 - 2f^ 



M 2 K± = mp^ + ^>: 2 a Q + v')zJ(Af K± -2f K0 -2f v 



M 2 K0 



^ (A' 2 + A 3 ^) Z m 2 (f K± + f K0 - 2f n ) (7f K0 - 5f K± - 2f w 
+ t\ (3A' 2 a + V) zj (4f I<0 - 2J K± - 2j n 



(5.2) 



12 

+ ^ (A' 2 + A 3 a 2 ) Z7 m 2 (J K± + J K o - 2f n ) (7J K± - 5/ K o - 2f n 



or 



and 



_ /2 1 ~ -^2 



12 

= ^Z m (W K± +W K0 + MU) (0.3) 

Ms = i(2M^ ± + 2Mio-M' ± )-iA;Z m 2 (7^ ± +7A--27. X2 
M 2 = (?n* + ^pWK -1 



M 2 = -V2 



2 

, Z P 



i (A4 ± + M^o - 2Ml ± ) (5.4) 



Vz" 1 (/ K± + f K0 - 2/J - ^A 2 Z" 2 (f K± + f K0 - 2// 



We will use (|5.2|), (|5.3|) to determine w , A 2 and v for given decay constants, M n ±, M K ± 



and Af^-o (once the wave function renormalizations are known). The parameters m' 2 and 

V' are independent of A3 whereas the dependence of A 2 on A3 is linear in A. Hence, the 
unflavored pseudoscalar mass matrix given by ( |5.4|) depends only to cubic or higher order 
on A3 and we will therefore neglect A3 in the pseudoscalar sector altogether. 

In the following we will make a first attempt to estimate the parameters Tn' 2 , V 1 and 

A 2 . We present the values to different orders in the quark masses in order to gain some 
intuition for the convergence of the quark mass expansion for these parameters. We use 
in this section the simplified kinetic term ( |4.9| ). To zeroth order in the quark masses the 



octet-singlet mixing vanishes and m 2 = 0. This yields the zeroth order relations for the 
mass of the rf and the scalar octet 

_i 3 Z m 

(5.5) 

OA 2 O q — + -M ,— . 



M 2 


3^_ 7 - X 


3 


Zm 

Z p 






2 

5 MJ 


Zp 
Z~h 
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Here we have used for the second identities renormalized couplings according to 
a = zia = \ + 7 K ±+ J '^o) -53.1 Me V 



v 



1 

Zm<y = 




_ 3 

Z m 2 v , 


2 y-l — 12 
m g — Zj m m g ) 


Z m 2 Ai , 


A2 = Z m 2 A 2 , - 



(5.6) 



Ai — Z m z \i , A 2 — Z m z A 2 , A 3 — Z m 3 X 3 
etc. This yields the zeroth order estimate 

z/°)— ~ 11500 MeV . (5.7) 
Z p 

A first possibility for an estimate of A2 uses the relation between the masses of the scalar 
and the pseudoscalar octets. For this purpose we need 



1 ,_^2 



- (2M 2 K , + M 2 a ) , (5.8) 

where Mr* — 1430 MeV. The ao-mesons, however, are not unambiguously identified. Usu- 
ally it is associated with the well established ao(980) resonance. On the other hand, the neigh- 
boring / (980) is often assumed to be an I = 0, KK bound state or a four quark state. One 
may therefore as well identify the ao(980) with the corresponding 1=1, KK "molecules". 
The only remaining candidate for the ao-mesons is then the possible ao(1320) resonance ||15|| . 
In this section we will use both possibilities with the notation M a ± = 1320(983) MeV. We 
obtain 

m h ~ 1394(1298) MeV 



x w_m _ 156.9(126.5) - 72.0 h/-l 

In addition to the uncertainty in the identification of the ao-meson the value of A2 depends 
rather sensitively on ratios of wave function renormalization constants. For this reason 
we will compute below the value of A2 from an expression involving only properties of the 
pseudoscalar mesons. 

To linear order in the quark masses the mass eigenvalues of the pseudoscalars are not 
affected by quark mass corrections to the kinetic term (see sect. §). We can therefore identify 
Mi = Mi and obtain 



1 (M 2 K± + M 2 R0 + Ml ± ) ~ (411.7 MeV) 2 (5.10) 



1 

12V = (jA 2 ao + ^) v = ^(M 2 K± +M 2 KQ -2Ml ± 

~ (388.9 MeV) 2 . (5.11) 
We may fix the parameter by the relation 

M 2 , = ml = 3 -^ua + 1^ (M 2 K± + M 2 R0 + m\ ± ) (5.12) 
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and find 



v (X)_m _ / 93 7 2 _ 2124 



MeV . 



The coupling A2 can now be infered from ( |5.11 



.(*) 



2 (M 2 K± + M 2 



K° 



2m! ± : 



3 (/** + f K o - 2f n ) a 
77.93 - 72.03 - 1 

V Zm 



2 Z v 2 (M^ ± + M^ + M'±; 

9 Z m al + 27 



(5.13) 



(5.14) 



We note that this determination of A2 uses the ratio of two quantities which are linear in A. 
The relation between f i and the meson decay constants fi is strongly influenced by quark 
mass corrections to the kinetic terms already to leading order in the quark mass expansion. 
From there we expect sizeable corrections to A2. In addition, an estimate of A2 requires 
information on Z p /Z m and involves differences in larger quantities. As a result, A2 will be 
poorly determined even once the quark mass corrections to the kinetic terms are included. 
We will find in section [13] typical values A2 — 15 — 30 even for Z p equal or somewhat smaller 
than one. 

Concerning the scalar meson masses to linear order in A we observe the relation 



nil 



Zm 2 

~z: m * 



z \ 1/2 
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(5.15) 



which translates with ( 5.11]) into 



2a (M 2 K± 



2 +M 2 



2M! ± 



ni 



Inserting experimental values this leads to the ratio 

-f- ^ 0.69(0.79) . 



+ 



m. 



(5.16) 



(5.17) 



These values are changed to Zh/Z m ~ 0.40 — 0.65 once quark mass corrections to the kinetic 
terms are included (see ( 11.6|) ). From the mass splitting within the scalar octet we infer 
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9o- A 2 + 3v , 



Ml* - M 2 



M Z K± + M 



2m!± 



0.67(2.38) . 



(5.18) 



This relation can be used to estimate the size of A3. 

IT,. 



Having computed the parameters u, m%, ctq, A2, A3, Z m /Zh, v and w from M^,, M^±, 



Mi , Mr±, Mi,, M ±, f%± and f n we can now derive other meson properties. Within the 

o d Q 

approximation Mj = Mj, f i = fi used in this section we discuss here briefly the non-flavored 
pseudoscalar mesons. For a determination of the pseudoscalar mixing angle 6 P to 0(A) and 
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the T) and rj' masses to 0(A 2 ) we need the off-diagonal element in the mass matrix for the 
neutral pseudoscalars to 0(A) 




+ f K0 -2f n ) (2\ 2 a 




(5.19) 



2/. 



- i fM 2 



+ M^o + M 2 ± 



We note that for Z p ~ J 
with Mg given by (|2.14|) and M 1 



the second term almost cancels the first one. Using ( |2 . 1 6| ) 
M v / ~ 957.8 MeV as an experimental input, one finds 
P = —18.7,-7.2,0.5 for Z p /Z m = 0.5,1.0,1.5. We will see in sect. ^| how quark mass 
corrections to the kinetic terms modify these relations substantially already to linear order 
in A. 



For the r\ mass we get to linear order the Gell-Mann-Okubo relation 



« 



Mi 



\{2M 



K± + 2M K0 - M n± ) ^ (566.3 MeV) 



(5.20) 



We can use ( [5.19|) together with the estimate of M g to quadratic order in A, ( |5.4|) , to 
compute a mass relation between M 2 and M 2 , to quadratic order in A. Inserting A2 from 
( [5. 14| ) and using M v / as an input we find 



521.2,535.7, 549.9 MeV 



(5.21) 



for Z p /Z m = 0.5,1.0,1.5. This is already very close to the experimental value M v = 
547.5 MeV as long as Z p /Z m is not too small. 

To summarize this section we find that the approximation of a quark mass independent 
kinetic term (|2TT|) or ( |4.9| ) gives already a reasonable overall picture of the scalar and pseu- 
doscalar mesons. The most important modifications from the quark mass corrections of the 
kinetic terms are expected for the value of A2 and the mixing angle 6 P . This will, in turn, 
influence the estimate of M 2 to second order in A and similarly the relation between M^, 
and Mp 3 . 



6 Quark mass corrections to kinetic terms 

The kinetic term ( |2.1| ) obtains corrections for nonvanishing quark masses. Expanding around 
$0 — diag(7J ) with (x. s ) = these corrections involve only the expectation value of h. We 
are only interested here in the kinetic terms for m and p. The most general corrections 
involving two derivatives and being linear in v and w can then be written in the form 

4£ = ^rn Tr hd^md^m + -u^p Tr h&*m . (6.1) 
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The term ~ u m leads to different wave function renormalization constants for pions, kaons 
and m 8 according to 

Z w = 1- uj m v 



(6.2) 



Z K ± = 1 + -u m (v + w) 

Z K o = 1 + -u m (v - w) 

Z s = 1 + uJ m v . 
This implies that the renormalized pion mass M w ± obeys 

Ml ± = Ml ± Z- 1 (6.3) 

2 

where M n ± is the mass computed in the previous sections, e.g. 

M 2 W± =m 2 m - l2 v + 6 7 (3v 2 + w 2 ) + 25 8 (v 2 + w 2 ) + 25 9 (v 2 - w 2 ) . (6.4) 

Similar relations hold for M|±, M 2 R0) Mf and Mf whereas Mf 8 = ~M 2 38 Z^Zp. There is 
also a mixed term ~ uo m wd^"mzd^m 8 . It gives corrections ~ w 2 to the 7r° and r\ masses and 
will be neglected here. 

One should note that the choice of Z n , Z K ±, Z K o and Z s is somewhat arbitrary. It 
depends on the convention for Z m , since a rescaling of Z m would result in a rescaling of Z^, 
Z K ±, Z K o and Z 8 . We employ here a convention for Z m where (cf. (|6.2|)) 



Z w + Z K ± + Z K o = 1 . (6.5) 
Neglecting corrections ~ w 2 we may then use 

1 / Zk±Zko \ 
UmV = 3 1 — Z 2 1 J ' ^ ^ 

The difference between M and M 2 influences the symmetry relations once expressed in 
terms of physical masses M 2 . In particular, we observe a modification of the relation fl3.4|) 
between Mf, M 2 ±, M R ± and M|- . Including the terms ( |6.1| ) this relation now reads (see 
eq. (|6.9| ) for the definition of /) 



Ml = ~ ( 2il4 ± %i + 2M| ^ - M 2 ± f^| - \\ 2 (J K± +J K0 - 2J? i . (6.7) 



2i 



8 



In consequence, the corrections due to the modification of the kinetic term influence the 
pseudoscalar mass eigenvalues to second order in the quark masses. Expanding fl6.7| ) to this 
order one finds neglecting terms ~ w 2 



Mo = 
8 3 



- [2M 2 K± + 2M 2 KQ - M 2 ±) - -u m v (M 2 K± + M 2 K0 - 2M 



l 2 (&* 



Uk± + ftco - 2/ w ) - -w m u (/ x ± + /jco + 4/ w 
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Here we have used that the deviation of Z %1 Zk±, etc. from unity also influences the 
relation between v, w and the decay constants f n , fx±, fx - The effect of the wave function 
renormalization on the meson decay constants is discussed in appendix [A] and leads to 

_j_ _i 

Z K ±fK± — Z R of K o = f K ± — f K o 

(6.9) 

i [z K lf K± + Z K lf K o - 2z;*f«) = \ (7a'± + 7x° - 2/J 

with 

f 7T = Z 7T 2 /tt 

and similarly for Jk±, Jk°- 

In summary, we will denote the physical meson masses and decay constants by Mj and 
/j. They correspond to a normalization of the fields with inverse propagator q 2 + M 2 in the 
vicinity of q 2 = —Mf. This is also the relevant field normalization for the decay constants 
- see appendix On the other hand, the quantities Mj and f i correspond to a common 
SU(3) symmetric wave function renormalization for the whole octet. Symmetry relations 
are therefore most easily expressed in terms of Mj and f\. For an approximation of the 
kinetic term to lowest order in the quark masses there is no difference between Mj and Mj 
or fi and f i . This approximation is sufficient to compute the meson masses to linear order 
in the quark masses, but not for decay constants and the mixing angle 9 P . In sections ||-|5| 
we employed a lowest order approximation to the kinetic term and therefore omitted the 
distinction between M and M or / and / for the quantitative estimates. On the other hand, 
the algebraic relations in the preceding sections are all expressed in terms of M and / and 
are therefore not altered by modifications of the kinetic terms. In consequence, the only 
necessary change for the quantitative estimates involves the relations between / and f ox M 
and M. 

We also note that the relations ( |6.2|) use a wave function renormalization which is defined 
for a common momentum q$ for the whole octet. In fact, one may view the derivative expan- 
sion as a Taylor expansion of the inverse propagators G~ 1 {q) around some fixed nonvanishing 
momentum q^ rather than an expansion around zero momentum, i.e. an expansion of the 
type G~ 1 (q) = M i + q 2 {Zi + 0(q 2 — q^)}- Here the depend on the choice of q$ through 
the normalization condition 

Zi = - 2 (G.\ql)-Gr\U)) 
M\ = G-\0). 

This condition, together with ( |6.5|) , also specifies the precise meaning of0 Z m . Similar def- 
initions also apply for Z p and Z^, but the momentum used can now be different from the 

8 For an expression of the propagators in terms of the fields <t> or 4> p the Zi stand for Z v ±Z m etc. and M i 

2 

should be replaced by M n ±Z m etc. 
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(6.10) 



(6.11) 



pseudoscalar octet momentum q$. For the definition of Z p it seems convenient to replace^] 
Qo by qf> = —^p, whereas for Zh one may use qf = —ml. We note that our definition fl6.11|) 
of Zi also specifies the precise meaning of the couplings u) m and u pm in ( |6.2j ). They multiply 
three-point functions at zero momentum for[^] h and momenta given by qf for p and m. All 
these specifications are irrelevant as long as only terms with two derivatives are included. 
The conceptual setting becomes crucial, however, once we go beyond this approximation and 
include higher derivative terms. 

In fact, the definition Mf = M i jZ{ yields the physical pole masses Mf only if we replace 
in ( |6.11| ) the common momentum qf by the individual locations of the poles at qf = —Mf 



This involves corrections of the inverse propagators ~ {qf — qf). We discuss these additional 
contributions from higher derivative terms more explicitly in section O. Here we note only 



the following general properties: Except for the mixing angle the higher derivative corrections 
can be absorbed completely in the wave function renormalization constants Z^. They lead to 
additional terms on the right hand side of ( |6.2| ) which are proportional to (qf — qf) according 
to 



GT\q) = M] + Zd + HM 2 ~ 9oV + \~Hf\q 2 - g 2 )V + 
= M] + Z l q> + H l {q 2 -q 2 l )q 2 + .... 



(6.12) 



Here the quantities Zi correspond to a normalization at q$ and are given by ( |6.2|) whereas 
the true wave function renormalizations Zi are defined at qf = —Mf, with 

Zi = Zi + AZ, 

1 12) ( 6 - 13 ) 
AZ t = H,{q 2 l -ql) + \H?\ql-ql? + ... . 

The difference qf — q^ is given by pseudoscalar mass differences and is therefore linear in 
the quark masses. More precisely, it is again linear in v and w. The Taylor expansion of 



G — G 1 (0) around q$ seems most reliable for y—qi somewhere inbetween the kaon and 
pion masses. We will choose the renormalization scale for Z m as 

= -~ (M 2 K± + M| + Mf ± ) (6.14) 

such that ( |6.5| ) remains valid. This implies that AZ n = —(AZ K ± + AZ K o) and we can 
absorb the higher derivative effects partially in a redefinition of an effective 

u m = u m (6.15) 



9 For convenience we will later also use = —Mf.. 

10 The normalization of h, however, is specified by Zh and therefore adapted to the behavior of the scalar 



propagator in the vicinity of its pole. 
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leading to 



Z K ± 



1 + 2 W m(« + w) 



(6.16) 



1 + -u m (v - w) 



Z K 

Z 8 = 1 + UJ m V + K 8 . 

In terms of uj m only Z 8 receives an additional correction 

K 8 = AZ 8 - 2AZ K . 

This yields the complete expression for M| to quadratic order in the quark masses 



(6.17) 



Mi 



~ (1 - K 8 ) (2M 2 K± + 2M 2 K0 - M 2 ± ) - hj m v (M 2 K± + 



M 2 



K° 



2Ml 



- J*. 



{Sk± + fxo - 2f w ) - -u m v (f K ± + f K o + Af n ) 



(6.18) 



For a computation of M v to quadratic order in the quark masses one needs, in addition, the 
complete expression for the mixing angle 9 P to linear order in the quark masses. 

In order to get a rough estimate of the size of the higher derivative corrections one 
may assume that the true inverse propagator does not deviate by more than 10% from the 
lowest order form M 2 + q 2 over a momentum range between — M 2 and —M 2 . This results 

in a typical bound H m ^ 0.1 \q 2 \ 1 and we expect the higher derivative corrections to be 
unimportant. Furthermore, within a systematic derivative expansion we can take to leading 
order a common H m for the whole octet. This yields 



M 2 K - Mi 



1 



:M1 



.3 A "3 

Inserting the leading order for M 2 on the right hand side of ( |6.19| ) we find 

K 8 ~ C(A 2 ) . 



(6.19) 



(6.20) 



We conclude that K 8 gives corrections to the pseudoscalar masses which are formally cubic 
in the m q . 

The deviation of Z^, Z K ±, etc. from one can lead to a sizeable change of the infered value 
for v and therefore to important modifications of the values of the couplings of the linear 
sigma model. Inserting ( |6.16| ) into ( |6.9| ) and ( |5.6| ) yields to leading order (neglecting isospin 
violation) 



Z m 



00 



Jk — S-k - 

Vk + u 



:U m v (f K + 2/ T ) 



U) 



(6.21) 
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We observe that the correction to the relation Q2.6Q is of the order (Zh/ 1 Z m ) l l 2 uj m (jQ. It is not 



suppressed by quark mass terms. On the other hand, the correction to a is of second order 
in the quark masses and therefore small. Taking uJ m v = —0.20 (see sect. 13) the relations 
( gglD result in 



7 



23.3 MeV, (Xo^ 53.8 MeV 



(6.22) 



to be compared with ( |2.7| ) and ([5.61). We conclude that the quark mass corrections to the 
kinetic terms are important for a quantitative understanding of the linear sigma model! 
These corrections also matter for a determination of the decay constant Jk°- We first notice 
the change in w induced by uo m v = —0.20 in ( 2.12 ). Using ( 6.1 6| ) we find 



w 



(fr) 



(Mf 



m 2 k± : 



3 

'2(l-uJ m v) {M 2 K -Ml) 



-0.67 MeV 



(6.23) 



This yields 



z h 



w 



1_ 



0.28 MeV (6.24) 



to be compared with (|2.12|) , /^o — /#± = 0.47 MeV, for ZJ m v = 0. The quark mass corrections 
to the kinetic terms reduce the isospin violation of the decay constants significantly! 

Let us next turn to the term ~ ^o pm d^p Tr hd^m in ( |6.1| ) which influences the octet-singlet 
mixing angle 9 P . This term leads to an off-diagonal kinetic term 



(6.25) 



Using m R8 = Z$ m 8 the inverse propagator for the fields p, is given in momentum space 
by the matrix 



G 



-i 



V 



z p (q 2 )q 2 + Ml 
^uj{q 2 )q 2 + Ml 



2 (&(q 



z 8 {q 2 )q 2 



2 V + M 2 p 



Mi 



(6.26) 



Here M| p is given to quadratic order in m q by fl5.4|), whereas M^ = m 2 , ( |5.12| ), and M| obeys 
( |6.7|) . We parameterize corrections from higher derivative terms for the diagonal elements 
of ( |6.26|) by the functions z p (q 2 ), z 8 (q 2 ). Our normalization of Z p and Z 8 is chosen such that 
z p (— M 2 ) = 1, z 8 (— M 2 ) = 1. Higher derivative corrections to the off-diagonal elements of 
( |6.26|) result in an effective g 2 -dependence of Co. For q 2 = —M 2 this correction is formally of 
third order in the quark masses, whereas for q 2 = —M 2 , it is of first order (since M^, — M^ is 
counted as being of zeroth order). In real life, however, M 2 , and M 2 are separated by a factor 
less than four. With the reasonable assumption that there is no dramatic g 2 -dependence 
of Cj we will treat uq 2 as a term linear in the quark masses for both — q 2 = M 2 and M 2 ,, 
and correspondingly count the higher derivative corrections as terms quadratic in the quark 
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masses. For a lowest order estimate they can therefore be neglected, z p (q 2 ) = z$(q 2 
Co(q 2 



LO. 



The diagonalization of ( |6.26|) can not be performed independently of q 2 anymore. As a 
consequence, the effective octet-singlet mixing angle 6L will depend on q 2 : 



tan 6 p {q z 



M 2 - M 2 + y/(Mj - M 2 ) 2 + 4(M| p + Coq 2 ) 2 Z, 1 

i 

2 



2(M 2 p + ooq 2 )Z 8 

The propagator in a diagonal basis corresponds to the eigenvalues of G, i.e. 



(6.27) 



G^W)=<l 2 + l{M 2 p +M 2 8 )±\ 1 -(M 1 



Ml 



M 



8 P , u 2 4 &M 2 p 2 



+ + — g 4 + 2 

Zj$ Zg Zj> 



The masses M 2 , and M 2 are given by the location of the poles of G for negative q 2 : 



(6.28) 



-1 r 



\(M 2 p +M i 



LO 



-Mi 



Hp 



1 



± {-(Mf-Mi) + 



M 



8p 



^M 2 8p (M 2 p + Ml) + 



a; 



(6.29) 



Expanding to quadratic order in the quark masses (with Co and M| p linear in m 9 ) one finds 
the relations 



M 2 + M 2 



M 2 , - M 2 



M 2 p + Ml - 2uMl p + 0? (M 2 p + Mg 
M 2 - M| + £ 2 (ilf 2 - Mg 



+ 



M 2 - Ml 



(6.30) 



Mj, - CoM 2 8p (M 2 + M 8 2 ) + uj 2 M 2 p Ml] (6.31) 



M 2 + M 2 - 2M 2 + 



2M| p 



M 2 + M 2 - 2M 2 



It is remarkable that the relation between M 2 , — M 2 and M 2 , + M 2 becomes independent of 
Co in this approximation. 

For an experimental determination of the octet-singlet mixing through the decay of rj' 
or r\ into two photons the relevant quantities will be p (q 2 = —M 2 ,) or 9 p (q 2 = —M 2 ), 
respectively. The mixing angle depends strongly on Co if CoM 2 , is of the same order of 
magnitude as Mf . This leads to a sizeable dependence on q 2 . If one intends to compute 
the octet-singlet mixing angle to quadratic order in the quark masses one needs in addition 
contributions to Co quadratic in m q . They arise from a modification of the kinetic term 
through 

4ii = l^' pm d,P Tr h 2 d»m + ... . (6.32) 
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This leads to a second order correction for Cj 



<Cj{q 2 ) = —\/3 UpmV + u. 



pm 



Uq 2 



(6.33) 



contains the higher derivative corrections where f w (—M 2 ) = 1 and / w (— ra^J 

from one has 
the formula ( 6.27|) remains 



where 

reflects the normalization of cj pm , u' pm at g 2 = — m 2 m . Also the deviation of Z- 
to be included to this order. We note that for z p (q 2 
valid if Cj is replaced by Cj(q 2 ), (|6.33|) . 



zs(q 2 



7 M v and My to quadratic order 



We are now ready to address the question of mass relations for the pseudoscalar mesons 
to quadratic order in the quark masses. First, we replace the coupling u, or equivalently 
M 2 , by M 2 , as a phenomenological input parameter. Our aim is to compute M 2 and the 



mixing angle 9 P as functions of M 2 ,, M 2 ± , M 



use the relations ( |6.31| ) and (|6.27|) which involve the quantities Mf , Mf p Z f 



K± , Mk , fx and /„-. For this purpose we 

-1/2 , * 7 -l/2 

and uZo 



The difference M 2 , 



M 2 is to quadratic order independent of Cj. The dependence on Z 8 



arises only indirectly through the correction ~ K$ in Mf (|6.18 ) and gives no correction to 
quadratic order in the quark masses, ( |6.20|) . The mixing angle depends already to linear 
order on Cj, but the difference between Z s and one is only needed to quadratic order. For 
the octet mass term Mf we use (|6.18 ) with 



U5 9 v 2 = -~\ 2 (f K± + f K0 - 2 f n ) 2 



(7.1) 



and determine A2 by (|5.14j) , with M 2 , replaced by the lowest order expression M 2 , + M 2 



I (2M 2 K± + 2M£o - M 2 ±). To the order relevant for our estimate M 2 p is given by 



Inserting this into the second expression ( |6.31|) one finally obtains the relation 



M 2 = \(2M\ 



x 



x 



r |± + 2M 2 K0 - M 



r U) + \ K ± + M| - 2M 2 ± ) 



■2/ w 




JK± + JK" + Jtt 



/x± + + fir 



9 (/*± + /j 



\ (3M 2 K± + 3M 2 K0 - 2M 2 ± Y 



(7.2) 



3 \Z, 



1 + 



- \ 2 

2M 2 :± + 2M 2 r0 - 1 r2 



(7, 



% - Ml 



(2Mj 



1± + 2M 2 K0 - M 2 , 



3M 2 - 2M| r± - 2M 2 + M 2 , 
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2 / Z m 

3 \ Z~ 



(M 2 K± + M| ) (2f K± + 2f K0 - f n ) - M 2 ± (f K± + f K : 



J (3M 2 - 2M 2 K± - 2M| + M^ ± ) (/ A - ± + ^ + / 
-K 8 (2M 2 K± + 2M 2 K o - M- 



— \ 2 



For cJ m = and K 8 = we recover the mass relation (|5.21| ) of sect. |5|. Based on ( |7.2| ) one 
can now evaluate the corrections ~ ZJ m and ~ K 8 quantitatively. Neglecting higher orders 
in {f K ± + f K o — 2fjj one has the approximate relations 

(7x±+7^Q-27,) ^ (/ ^ ± + f K o - 2/0 _ 1_ ^ 



+ /*<> + / J + fK° + /-) 
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In the following we also drop K 8 since it is formally of third order in the quark masses. If 
we linearize (|7T2| ) in ZJ m v = § 1 u m (f K ± + f K o - 2f n ^ and - l) we find 



M v = (535.7 MeV) 




0.511 u m v 



(7.4) 



Within the range Z p /Z m = 1.0 ± 0.2 we see that the corrections ~ (Z p /Z m — 1) only amount 
to at most 7 MeV. For a rough estimate of the size of UJ m v we can neglect them. Comparison 
of (|7.4|) with the measured r\ mass yields for Z p = Z m 

u m v ~ -0.043 . (7.5) 

We should note, however, that the linearization in uJ m v is not reliable anymore for uJ m v ^ 
—0.15. This is demonstrated in fig. [I] where we plot as a function of u m v for various 
values of Z p /Z m . For this plot we have evaluated all matrix elements in ( |6.26 ) to second 



1/2 

order in A (keeping the full cU m f-dependence, e.g., in Z 8 , though) and diagonalized the 
matrix without further expansion in A. We have neglected higher derivative contributions 
which can not be absorbed into the wave function renormalizations, i.e. we used z p = z 8 = 1, 
f u = 1. For uj we have used ( |8.14| ) and u m = uj m . The nonlinear effects due to terms 



which are formally ~ 0(A 3 ) and higher are reflected by the deviation of the curves in fig. [1] 
from the tangents at u m v = 0. Because of the important nonlinearities for UJ m v ^ —0.15 
one finds that there is a second solution with Z p ~ Z m , namely for u m v ~ —0.22. For 
large values of \u m v\ the rj-rf mixing starts playing an important role. We will see below 
that realistic values for the decay rates rj —>■ 2 / y and v( — > 27 are only consistent with 
this second solution for uj m v. In this region the formal quark mass expansion does not 
converge well anymore! This breakdown of the quark mass expansion is even more apparent 
in the relation between M 2 , and M 2 since linear and higher order mixing effects occur 
with the same sign for this quantity. We plot in figure |2| the value of M p as a function 
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III 

500 1 ' ' I 

-0.25 -0.2 -0.15 -0.1 -0.05 0.05 0.1 



UJ m V 



Figure 1: The plot shows M n as a function of tv m v for various values of Z p /Z m and u m — ZJ m . The solid 
line corresponds to Z p /Z m = 1 and the difference in Z p /Z m between two adjacent lines is 0.1. The horizontal 
dotted line indicates the experimental value M n ~ 547.5 McV. 




-0^25 -0.2 ' -0.15 -o'.l -o!o5 0.05 ' 0.1 



Figure 2: The plot shows M p as a function of uj m v for given M,y = 957.8 MeV, uj m = uj m and Z p /Z, 
varying between 0.7 and 1.3 in steps of 0.1. The solid line corresponds to Z p /Z m = 1.0. 
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of uj m v for a fixed value of M v > = 957.8 MeV and various values of Z p /Z m neglecting all 
higher derivative corrections to the kinetic terms. As a first observation one sees that the 
dependence of M p on Z p /Z m is rather week. Furthermore, for u m v ~ —0.22 the ratio M 2 /M 2 , 
has decreased to about 0.77 despite the fact that this effect is formally of second order in A 
if M 2 is counted as (9(1). A partial explanation of this strong mixing effect is related to the 
observation that for u m v ~ —0.22 the value M p ~ 839 MeV is actually almost comparable 
to M 8 = 579 MeV. (The values are for Z p /Z m = 1.0.) A counting where Mg/M 2 = 
becomes therefore doubtful. As an alternative, we may count the anomaly contribution to 
Mp, i.e. MpZp/Zm — rr? m ~ (742 MeV) 2 also as an effect of order A since its size is not too 
different from the £77(3) breaking induced by the mass of the strange quark. In this counting 
scheme all elements of the rj-rj' mass matrix have the same order of magnitude. Eigenvalues 
to a given order in A involve then all matrix elements to this order and the mixing effects 
are naturally large. Nevertheless, since M 8 /M 2 = 0.69 is still substantially smaller than one 
the real situation is somewhere in the transition region between the two counting rules. Our 
approximations to order A 2 are consistent with both ways of counting. 

In the remaining sections of this paper we will gradually collect information on the 
quantities Z m /Z p , uJ m v, u and K 8 . We should mention that the quantities uJ m v and K 8 only 
involve properties of the effective action for the pseudoscalar octet. There is, in principle, 
no information which goes beyond the one contained in chiral perturbation theory. The 
quantities Z m /Z p and u> also involve the rf and go beyond standard chiral perturbation 
theory. Of course, if one were able to predict directly quantities like ZJ m v one would gain in 
addition information on some parameters appearing in chiral perturbation theory. 



8 Kinetic terms in the linear meson model 

In this section we discuss in more detail the derivative terms in the context of the linear a— 
model. Our first aim is to gain information about the size of Z m /Z p , uj m and u. Expanding 
( |4.6| ) in powers of $ p , $ s and Xp we observe that only those terms contribute which have 
both derivatives acting on pseudoscalar fields. Also all fields without derivatives must be 
scalars. We can therefore replace in 



i i (8 "^ 



and for the fields without derivatives 



$ — ► a + -=$ s 

f (8.2) 



The terms appearing in £ kin ( |4.6j ) have been selected such that the leading terms contributing 



to Z p /Z m} uj m and uj pm have been included. For this purpose we first classify the invariants 
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contributing to terms with two derivatives acting on pseudoscalars which are at most linear 
in $ s . The reasoning is somewhat lengthy but straightforward: By partial integration all 
invariants involve at most one derivative acting on a given <3>. We can then distribute the 
two derivatives either within the same index contraction (with 5 ab or e abc ) with respect to 
SUl(3) x SUr(3) or among two different such structures. Traces involving six and more 
powers of $ have at least one combination ( or $$t) without derivatives in the chain. 
By a suitable definition of the invariants this may be replaced by — | Tr($^$)^ and 

such invariants contribute therefore only to higher order in the quark masses. For instance, 
combinations involving two factors of or $$t (without derivatives) contribute at most 
to order <£>g. The only invariant involving a trace of six factors $^ or $ that may contribute 
linearly in $ s is therefore the term ~ W v . There are also two terms quartic in $ with 
couplings X£, X~. To linear order in <3> s the structures involving e-tensors are the one ~ U v 
and terms not listed in (|4.6|), namely 



+ T^d^ aibl ^ a2b2 (d^bMta^atbs + ®a 3 bM ia A®a4b 3 ~ ^$0363 Tr 

+ T^d^ aibl $ a2b2 U a3b4 d^l 4a4 $ a4b3 + ld^ asb3 Tr (8.3) 



3 

+ TW U aibl ^ a2b2 d^ a , b M a d^ aibz + U aibl ^ a2b2 ^ asb3 




J pm- 



3 

* ± ip ^0161^0262 \y ^a3D4 L/ M^6 4 a4^a-403 ~ ^0364°' ^6404^^0463 



T "^0364 ^6404 °>^ a 4f>3 



These invariants do not contribute to 

or but they may contribute to to m or Up 
We treat C^iT) as a higher order correction to the term ~ U v and neglect this piece in 
the following. Next we turn to the case where the two derivatives act within two different 
SUl(3) x SUr(3) invariant index structures. Since the invariants p, r 2 , r 3 and £ are at least 
quadratic in the pseudoscalar fields the terms ~ d^pd^p, d^r 2 d' 1 T2, d^d^r^, or d^d^ do 
not contribute to Z m , Z p , uj m or u>. The same holds for mixed terms like d^pd^^ and for 
index structures of the type Tr $t<9 M $ Tr etc. On the other hand, the pseudoscalar 

invariant uo = i(det $ — det $') contains a term linear in the pseudoscalar fields 

d^u = -V6a 2 d^ Xp + Tr $% X p + Tr $ s d^ p - i= Tr $^$ p + . . . . (8.4) 

More generally, we can construct invariants with a possible contribution to Z p /Z m , u m or 
by Lorentz contraction of CP-odd factors like d^u, 

^ e aia 2 a 3 bib 2 b 3 $ , $t <j, _ » <T>t <j,t $t \ 

tc ' ^^^0161^0262^0364^6404^1463 ty /i^ r ai6i ^0262 0364 ^0404 ^0463/ ' 

ie aio 2 03 e 6i626 3 ^ {$ ai6l $ a262 $ a3b4 $t 4a4 $ a463 - $ 0l 6 1 <62<64^64a4<6 3 } , 

4 e aia 2 a 3 b 1 b 2 b 3 $t $ i - $ f , $ f , $ f , 9 $i $t , 1 

tC ' |^aiOl ^Q202 ^1364^^6404 a 4°3 ^0161^0262 0364 A» 6 4 a 4 0463 J 
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Within these structures we may replace subsets $$t<fj by the combination — 1$ Tr 

which is ~ $ s . We will keep here only the term which is ~ d^ud^uu and consider the other 
contractions as higher order corrections to the term ~ V v . The phenomenological analysis 
below will indicate that the "anomaly" terms ~ U v and ~ V v are not very important. This 
justifies to neglect corrections to them. We have also included in ( [4.6| ) some invariants which 
do not contribute to Z p /Z m , uj m and u> pm . This is partly for the purpose to demonstrate 
that not much can be learned about the ratio Zh/Z m from exploiting the symmetries of the 
linear sigma model. 

The contribution of the term ~ V v to Z m , Z p and to can be read off from 

+ s/3V v a 2 d» X pTr + (8.5) 

Next we turn to the terms ~ X~ and in ( |4.(j| ). Whereas gives no contribution linear 
in the quark masses the term ~ X~ yields 



Adn(X;) = ^X-^ {Tr^ p d^ p + d^ X pd^Xp) + ^X-a Tr^%d^ t 
+ l=X-a d^ Xp Tr $ s d,% + -Lx-d^x P Tr . 



■6) 



Similarly, one obtains 

£kin(Wg = ^Wy Tr$ s d»%d^ + ^W^d» X pTr$sd^ p 

f 3 1 
+ W 9 o\ l-Tr^d^pd^p- —Tr^ 2 s Trd^%d^ p (8.7) 
^ 8 24 

+ I Tr $ a d»$ p $ a d^ p + Tr + 1 Tr $^xAx P } + • • • • 

For the piece ~ U v we use 

c t a\b\" a^bi® 0363 ~T~ "ai 62 ^12^3 "^3^1 "t" "aif>3 ( Ja2&i'Ja3&2 



and find 



® a\b\® a,2b?P a-sbi "0163^0262 0361 "0162 &2&1 03&3 



- Tr $ s 0"$ p fy$ p + ^rf^Xp Tr $ s ^$ p . 



•9) 



In our truncation no other second derivative invariant (except the one ~ Z v ) contributes 
to Z m , Z p uj m , uJ pm or uJ pm - Combining ( |8.5| ), ( |8.6| ) and ( |3.9| ) with the term ~ Z^ and 
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comparing with ( f4.9|) , (|Q|), ( |6.32|) and ( |6.33| ) yields 

Zm Z^ 



(Or. 



10. 



pin 



10, 



pm 



10 



+ X'al + QV^ - 2U^a 
X~a + l -W^l-U^)Z h ^Z rt 

1 



X-a + ^=W v o* - VdVyd* + ±=u)j Z p *Z h *zJ 



w^i + ^-v^\z P *z h l z 



f K± + f — ^fjr 

J_ K± + Jko - 2J W (J_ x 5 
f K ±+f K °+L bV2 * ° 6^2 



^ + ^X^ + ^Wy -V2V^ 



Zn 2 Z„ 



(8.10) 



We observe the relation 



UJpm 



1 
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\Z P J cr \Z } 



z p \- 



(8.11) 



which is independent of the values of X v , W v , U v and V v . It leads to an estimate of Cj to 
lowest order in the quark masses: 



U! 



(1) 



M 2 ^ 



2a 



1 - w m v 



(8.12) 



If we furthermore assume |6V^ctq| , 1 1 W^cr^j <C X (see section [10]) we findQ (for V v — 0, 



1 



to, 



7 \ 2 7 



pm 



2\^6a \Z P J Z, 



h 



0J r . 



This leads to an estimate to second order in A 

1 



Co® = w« 



2h 



2"\/2 I 2p 



z„ \* + J_f 1 _z ]L 

3o"o V Z m 



6(Jq \ Zj m 



(8.13) 



(8.14) 



These results can be used for a computation of the rj-rf mixing angle to second order in 
the quark masses. If we neglect the higher derivative corrections (p m = u m , K 8 = 0, 
z p — z 8 — fui = 1) the mixing angles 6 p {rj) and 9 p (rj') relevant for the two photon decay of 
the rj and rj', respectively, depend on two parameters ui m v and Z p /Z m — 1. We plot in fig. |] 
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Figure 3: The plots show the pseudoscalar mixing angles 9 p (rf) and p {r}') as functions of uj m v for various 
values of Z p /Z m and u> m — oJ m . The solid line corresponds to Z p /Z m = 1 and the difference in Z p /Z m 
between two adjacent lines is 0.1. 



these quantities as functions of u m v for various values of Z v jZ m — 1. Assuming 



Zm 



LO m V 



2(/k-A) 



and comparing this value with 

we find that the first contribution is smaller than the second one. In the approximation 
Z p = Z m the quantity Cj is simply related to oo m v by 



LO 



(1) 



LOmV 



LO 



(2) 



LO m V 



1 + 



z h 



4a 



5.15) 



The existence of a relation between Cj and u> m is crucial for the predictive power of our model 
since it is necessary in order to relate the rf-rf mixing angle to other observables. Within our 
truncation and to lowest order in A the relation (|8.12|) is a pure symmetry relation without 



any assumption on the values of the couplings , W v , U v , V v \ In contrast, the kinetic 
term of the scalar mesons is independent of V v , or X~, whereas Z h receives contributions 
~ etc. No new relations can be obtained in this way. 

We observe that the couplings U v and V v violate the axial C/a(1) symmetry and are 
therefore connected to effects from the axial anomaly. Anomaly contributions to kinetic 
terms of the pseudoscalar octet are often counted as higher order corrections. If we decide 
to do so we obtain the leading order relations 



(8.16) 



In this approximation only the kinetic term ~ X^ is essential. The importance of X 
is also manifest in the leading mixing approximation. Only this term is generated by the 

and Z m — 



11 The coupling uj' pm cannot be expressed in terms of 
invariants are linearly independent. 



This reflects that the kinetic 
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mixing with the divergence of the axialvector field <9 M p^ ("partial Higgs effect"). In fact, 
the effective action for the (pseudo)scalars contained in $ receives contributions from the 
exchange of other particles. Prominent candidates are, of course, the vector and axialvector 
fields. We have computed in appendix |B| the contributions from the exchange of the vector 
and axialvector meson octets as well as the associated h and + ~ states corresponding to 
the divergence of the (axial)vector fields. Up to order <3> 4 this exchange only contributes to 
derivative terms. In terms of the couplings which contribute to the kinetic terms to linear 
order in A only X~ gets a negative contribution from the exchange of the h state whereas 
all other couplings remain unaffected. We find a large effect 

u^v ~ -0.15 . (8.17) 

Comparison of this value with figs. [TJ-Q indicates a large 77-77' mixing in a range where 
nonlinear effects in u m v are already important! 



9 Decay constants of 77 and rj 



The decay constants f n o, f v and f v > are experimentally determined from the partial decay 
width of the tt°, rj and 7/ into two photons 

a 2 M 3 

r(,-2 7) = ^ (9.1) 

V 



and similarly for 7/ and tt°. The experimental values for the decay widths are [pTS 

r(7r°^ 2 7 ) = (7.78 ± 0.56) eV 

T(77^ 2 7 ) = (0.46 ± 0.04) keV 

T(rj' ^ 2 7 ) = (4.26 ± 0.19) keV (9.2) 
which yields (/ w = (92.4 ± 0.3) MeV) 



(fv 

u 



exp 



cxp 



cxp 



1.00 ±0.04 
1.06 ±0.05 

0.81 ±0.02. (9.3) 



To lowest order in the quark mass expansion the decay constant^ f v is related to fx = fw 
by SU(3) symmetry (cf. appendix 0) 

4 0) = V3f v8 = V3U . (9.4) 



12 We use for f n and /,,' the conventions of JL9[ , with the warning that factors y/S and -y/3/8 appear here 
as compared to a perhaps more natural convention based on SU(3) symmetry. In the limit of exact SU(3) 
symmetry the quantity corresponding to is J,^. 
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Similarly, for Z v = Z m one finds for f v > 



Jrj' 



(0) 




(9.5) 



The experimental values Q9.3| ) differ substantially from this estimate. In appendix [A| we 
have computed (( |A.24j )-( |A.26| )) the corrections to f v8 and f v > as well as the corresponding 
constants for the singlet current in the r\ (f v o) and the octet current in the rj' (f v 's)- Ex- 
panding the Z 1 ! 2 factors to linear order in ZJ m v and neglecting additional higher derivative 
corrections one obtains the ratios 



fyS 
fir 



fy'O 



1 



2f K ± + 2f K o - f n + ]-u m v (f K ± + f K o - 2f n ) 



+ ~K S {2f K ± + 2f K o - f v ) 



(9.6) 



1 

3^ 



fx± + Ik" + fn - -u m v Uk± + fxo - 2f n ) 



Z v \- 



Extracting the terms linear in the quark masses and neglecting isospin violation yields 

4/*: - U 



frfi 

u 



3U 



1.3 



f 



r/'O 



2f K + U (Z£$~ x . Xh (h- 



3/. 



(9.7) 
(9.8) 



We note that this estimate for /^ 8 / f n agrees well with the one obtained in chiral perturbation 
theory (U/U) x pt ^ 1-25 % p. Denoting by 6 p ( V ) = 6 p (q 2 = -M 2 ) and 9 p (rf) = 6 p (q 2 = 
—M 2 i) (see (|6.27|) ) the octet-singlet mixing angles relevant for the two photon decay of the 
i] and r/', respectively, one can compute pC§ the effective decay constants for these decays 
from 



1 1 fcos6 p (ri) y/8sm8 p (rj) 

fn V% \ fvs 

1_ = 1 f sm9 p ( V ') | 

fn' \ frf 8 




(9.9) 
(9.10) 



Here we have neglected the mixing of ttq with r\ and r/' which induces corrections to the decay 
constants ~ w 2 . Since also Z n o equals Z n up to corrections ~ w 2 we will not distinguish 
between f n o and f n in this paper. 

We plot in fig. £| the decay constants f v and /,,/ as functions of u m v for various values 
of Z p /Z m , assuming uJ m = u m . It is obvious from these plots that no satisfactory solution 
exists for the value uo m v ~ —0.05, corresponding to the first solution for M n (cf. fig. P which 
remains within the range of convergence of the quark mass expansion. On the other hand, 
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Figure 4: The plots show the ratios f v / f n and fn'/f-n as functions of Lo m v for various values of Z p /Z m and 
u) m = ZJ m . The solid lines correspond to Z p /Z m = 1 and the difference in Z p /Z m between two adjacent lines 
is 0.1. The experimentally allowed windows (la) for both quantities are bounded by the horizontal solid 
lines. 
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it is encouraging that for the second solution u m v ~ —0.22 [Z p jZ m = 1) f v is already very 
close to the experimentally allowed window! The deviation of /,,/ from its experimental value 
by around 25% is not completely unexpected for the approximations employed so far. In 
fact, already the uncertainty of a first order computation of f v and f v i in powers of quark 
masses should be of the order of 10%. On top of this, the neglect of higher derivative terms is 
less accurate for q 2 = —M 2 , and one expects a less convergent expansion for f v >. If Z p /Z m is 
treated as a second free parameter we see that u m v = —0.20, Z p /Z m = 0.9 provides a solution 
for which f v and f v > are within 10% of the experimentally allowed windows. Furthermore, 
these values of u m v are quite close to the ones estimated in appendix [B] from the exchange 
of higher _+ states (cf. QS.17|) )! We conclude that all observations fit together in a picture 
with large mixing in the 77-77' sector. As discussed at the end of section [7] the anomaly term 
comes out relatively small in this case. The naive quark mass expansion is then expected to 
converge well only for the flavored mesons whereas its convergence is unsatisfactory in the 
77—77' sector. There it can be replaced by a modified expansion where M p is also counted as 
O(A). 

Apparently, the kinetic term ~ which is induced by the mixing with the higher + 
multiplet plays a very important role in our picture. One is tempted to assume that this 
term dominates the rich structure in the kinetic terms for the pseudoscalars. The hypothesis 
that all deviations from a standard kinetic term for the pseudoscalars (i.e., eq. (|2.1|) with 
Z p = Z m ) are due to mixing leads to a highly predictive scheme. In the limit where X~ is 
independent of momentum this leads to Z p ~ Z m , ZJ m ~ tu m . It is far from trivial that for 
Z p ~ Z m and UJ m ~ u m there exists a value uj m v = —0.20 for which the quantities M v , f v 
and f v ' are all compatible with observation! The "robustness" of these "predictions" can be 
estimated from table || in sect. [13] where we also give numbers for various values of Z p /Z m 
and uj m /u) m . A more detailed study of the effect of mixing with the higher h states can 
be found in appendix 0. In particular, the coupling X~ becomes effectively momentum 
dependent due to propagator effects. This leads to 

(i) a contribution to a g 4 -kinetic term and therefore to (ZJ m — tu m )/uj m ~ 0.1 (cf. ( |12.3|) ); 

(ii) an effective momentum dependence of Cj which gets multiplied by the factor f w = 

{M 2 P - Ml) I [Ml + q 2 ) with M P > 2000 MeV; 



(iii) similar contributions to z p and z 8 , ( |C.13|) in ( 6.26 ); 

(iv) a correction to ( |A.26| ) given by ( |C.14] ); 

(v) an effective Z p (normalized here for q 2 = —M 2 ) obeying Z p /Z rn = 0.99. 

Due to the large uncertainty in the value of Mp relevant in the 77 — 77' sector we have not 
included these higher derivative effects in the figs. [l|-f|. A quantitative discussion can be 
found in sections [L2| (fig. ||) , [L3] and appendix 0. 
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10 Expansion in the chiral condensate 



There are various scales characteristic for the amount of spontaneous chiral symmetry break- 
ing: the chiral condensate | (ipip) \ — 200 MeV, the pion decay constant f n ~ 90 MeV and 
the constituent quark mass m q ~ 300 MeV. All these scales are typically smaller than the 
characteristic scale for the formation of the mesonic bound states, k v ^ 600 MeV ||, |6| or 
typical mesonic mass scales unrelated to the Goldstone phenomenon — the latter being 
around 1 GeV with the lowest one given 

by0 M p o = 770 MeV. The question emerges if the 
typical scale appearing in the parameters 77, X~ , A 3 , etc. is larger than characteristic scales 
associated with chiral symmetry breaking. Could it be that chiral symmetry breaking is 
related to an additional small parameter leading to a suppression of contributions with high 
powers of <To? The existence of a small parameter associated to cr would enhance the predic- 
tive power of the linear sigma model since the terms with lower powers of do would dominate. 
Together with a derivative expansion it would allow to classify invariants according to their 
dimension with a suppression of higher dimensional operators. 

Within the linear sigma model we observe distinct mass scales of very different origin: 
Whereas the expectation value o"o measures the strength of spontaneous chiral symmetry 
breaking, the scale v indicates the size of the explicit breaking of the axial Ua(1) because 
of the chiral anomaly. Furthermore, there are hadronic mass scales which are not directly 
related to chiral symmetry breaking or the axial anomaly as, for example, the string tension 
or the glueball masses. One expects that the last type of scales dominates the dimensionfull 
parameters in the effective potential and the kinetic terms in the limit of vanishing anomaly 
and o"o = 0. Since for mass terms etc. <jq is multiplied by some dimensionless coupling 
constant, a typical parameter for testing the convergence of a cr -expansion could be x a = 
^2<Jq/v. Using the values of the parameters determined in section^ yields x a ~ 0.3 whereas 
including the quark mass corrections to the kinetic terms gives x a ~ 0.2. This seems indeed 
to allow for the possibility that an expansion in powers of <To does not converge too badly. 
We will see that this picture is confirmed by the size of other dimensionless ratios involving 
powers of oq. 

In order to make a guess for the size of contributions with higher powers of Oq we first 
note that the smallness of cr is partly due to the smallness of Z m ||. Since the physics 
cannot depend on the choice of the scaling for the field only ratios which are independent 
of Z m can appear in measurable quantities. This includes combinations like X^^jv = 
\ 2 (?o/v, A 3 ajj/z7 = A 3 cx^/V, X~al/Z m , U^q/Z^ V v a^/Z m . The physical scales hidden in 
the renormalized parameters can be better appreciated if we choose (somewhat arbitrarily) 
a fixed Zffl such that the dimensionless couplings are of order one, say Z<£> = 0.15 such 

that A^ ~ A 2 /50. With this scaling of the field one has cif } = <ro(Z$)~ 1/2 - 137 MeV 
and = u{Z^Y^ 2 ~ 540 MeV. We will now assume that dimensionfull parameters 

like A 3 ° are given by powers of a characteristic scale which we take to be around 700 MeV, 

i.e. A~3 0) ~ (700 MeV)" 2 , ~ (700MeV)~ 2 . This suggests typical ratios |A 3 ^/z/| ~ 0.015, 

|A 3( x 2 /A 2 | ~ 0.08, \X~a 2 /Z m \ ~ 0.25, \U v a /Z m \ ~ 1.3, \V v a%/Z m \ ~ 0.01. Obviously, these 

13 The p-meson is perhaps somewhat special because of approximate gauge symmetry, see appendix H. 
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numbers can only be used as rough guesses. There may be additional small coefficients - 
this is obviously necessary for the contribution ~ U v if Z p is in the vicinity of Z m (see flS.lOD 



- or relatively large group theoretical or dynamical factors. If the contribution from 
does not dominate u m we may estimate \X~a1/ Z m \ ~ |a; m Z^ Oo| — 0.5 (for u m v = —0.20) 
which is somewhat larger but still compatible with the above guess. We conclude that the 
a"o-expansion converges at best slowly. For low powers of <To group theoretical factors or 
dynamically small quantities (i.e. Z$, U^Wq) remain very relevant. Nevertheless, we find 
it very unlikely that terms with high powers of gq dominate those with low powers. Even 
the conservative assumption that terms with high powers of <To are bounded in size by the 
strength of terms with lower powers has important implications! 

As an example we compare the contributions ~ §00-^2 and ~ 3<XqA3 to the cubic coupling 



76 Q4.34 ) which determines the mass split within the scalar octet. If we assume that the 



second term does not exceed in size the first one we obtain the bound 

7e < \ (18a A 2 - v) (???) 2 ~ 17.1 GeV (10.1) 
where we have used A 2 = 21.3, v = 6447 MeV and Z h /Z m = 0.35 (see sect. [13]). With the 



help of ( [2-19| ) this can be transformed into a bound for the mass difference between the Kq 
and do mesons in the scalar octet 

M 2 Kg - M\ o = 3-f 6 v < 0.7 GeV 2 . (10.2) 

This relatively conservative bound seems to disfavor the interpretation of the ao(980) reso- 
nance as a member of the same octet as the Kq(1430), since in this case the difference in 
mass squared would have to exceed 1 GeV 2 . In simpler words, it seems at first sight unlikely 
that a strange quark mass of about 180 MeV produces a mass difference between strange 
and non-strange scalar mesons of 450 MeV. Yet, we notice that ( |10.2| ) is subject to quark 



mass corrections from kinetic terms, and we will come back to this issue in sect. pLlj . 

Before closing this section let us comment on the question if the limiting case — > 
can be used as an expansion point within a generalized class of linear sigma models. (Of 
course, the sigma model corresponding to low-energy QCD has a fixed value of cr Q .) Let 
us consider an effective quark-meson theory which is supposed to be valid for momentum 
scales below some cutoff (or compositeness) scale k v . The "classical action" of such a model is 
parameterized by a potential and, in particular, a mass term W^ik^). (Quantum fluctuations 
of modes with momenta q 2 < fc 2 change the form of the effective action and lead to an 
effective potential as parameterized by ( |4.22| ) ) . The size of ctq can be influenced by the 
meson mass term fn 2 (k ip ) at the scale k v . If the phase transition associated to a variation of 
m 2 (/C(p) were of second order the order parameter <To could be arbitrarily small. An expansion 
in powers of <Jq would then always be meaningful for small enough^ ctq. For three flavors 

14 At nonvanishing temperature <Jo(T) typically decreases as T increases. For a second order high tem- 
perature phase transition one can always expand such a system for a very small even though perhaps not 
arbitrarily small value of a®. In the limit Co — > one may encounter nonanalyticities associated to the 
critical three dimensional behavior at the transition. 
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the anomaly induces a first order transition and excludes arbitrarily small values of a®. 
Nevertheless, for small ctq a polynomial expansion of U (a) should be meaningful and we may 



stop after the term quartic in a. In the limit of equal quark masses the potential ( 4.22 ) gives 



2' 



U = Sm 2 g a 2 - -vol + -A lC r 4 + ( 3m] + -va - 9X^1 ) a 2 - vo z + -Aicx 4 . (10.3) 



The requirement U(a ) < U(0) implies a lower bound 



9Ai 

whereas the positivity of the mass term at o"o requires 



^ 2 >7rr (^o + emf) (io.4) 



°o > ^ H - 2mJJ . (10.5) 

On the other hand, the dimensionless coupling Ai is typically bounded from above as a result 
of the "triviality" of <3> 4 -theory. (More precisely, the infrared interval of allowed renormalized 
quartic couplings is bounded.) Comparing ( |10.5|) with the definition of x a we find that for 
vanishing quark masses {m 2 = 0) the expansion coefficient must obey x a > and can 
therefore not be arbitrarily small. Despite this caveat there seems to be enough room for a 
meaningful a"o-expansion. It is interesting to note that for given o"o the inequalities ( |10.4[ ) 



and ( |10.5| ) can also be used to establish lower bounds for Ai which hold for a polynomial 
approximation ( 10. 3j ). Taking z/, m 2 and o"o from section ( |T3"D one finds Ai > 48.9 and 



Ai > 1.1, respectively. 



11 Mass relations for the scalar octet 

In this section we want to determine the masses of the members of the ++ octet contained 
in $. Together with the ++ singlet these scalars play for spontaneous chiral symmetry 
breaking the same role as the Higgs scalar in the electroweak theory. Since the isospin 
singlet member of the octet has the same quantum numbers as a possible scalar glueball the 
determination of its mass is also important for the identification of glueball state candidates. 
We restrict most of the discussion to the linear order in an expansion in powers of quark 
masses and we neglect isospin violation. The quark mass corrections to the kinetic terms for 
the scalar octet to linear order in A arise from an interaction analogous to (|6.1|) 



41; S) = ±u h Trh& i hdph. (11.1) 



Since in the linear sigma model several of the generalized kinetic terms (KM) contribute to 



this invariant we will treat uj^v as a free parameter to order A. In addition, we will consider 
a particular term contributing to second order in A 

& = 2ChTr{(hd»h-& t hh)(hd ll h-d ll hh)} . (11.2) 
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This term is induced by the exchange of the scalar state contained in the divergence of the 
vector meson field d^py (cf. appendix |B]) with a sizeable coefficient (h^ 2 ^ 0.02. Combin- 
ing ( 11.1 ) and ( 11.2 ) and neglecting terms ~ u^hV 3 this leads to different wave function 
renormalizations for a , K* and /g 



Za — 

Z K * = 
Z fs = 1 + u h v . 



uj h v 



1 + -u h v ) ( 1 - % h v' 



;ii.3) 



We next have to include the effects of higher derivative terms using the definition of the wave 
function renormalization constants ( |6.11|) . The discussion is completely analogous to sect. |6] 
and the dominant higher derivative terms lead to a replacement of ujh by ZJh = c^+2A.Z#*/i>. 
(Here we have adopted a definition of Zh such that 2Zk* + Z ao = 1 for (h = 0.) We note 
that the effects from a nonvanishing Ch can be absorbed in an effective mass term 



M 



i<:: 



M 



;ii.4) 



Up to the replacement of the physical mass M\, by the (h dependent quantity M\* 
discussion systematically only includes terms linearQ in A. 
From (15.81) we now obtain 



our 



™l 



2M K , + M n 



1 r 



2M£. + Mi o + ZJ h v [Mk* ~ Ml 



;il.5) 



Hence, the linear quark mass corrections to the kinetic terms ~ Uh modify ml only to 
quadratic order in A. For M ao = 1320(983) MeV, UJ h = and ( h = we findF] m h ~ 
1394(1298) MeV. The dominant correction to the lowest order relation is most likely due 
to the term ~ ( h . For ( h v 2 = 0.02 one obtains m h = 1303(1200) MeV whereas ( h v 2 = 0.04 
yields already a large shift to rrih = 1206(1093) MeV. We are now in a position to compute 
Zh/Z m from ( |5.16| ). For UJ m v = —0.20 we obtain 



0.39(0.45) , 0.45(0.53) , 0.53(0.64) 



(11.6) 



where the three values correspond to (h v2 = 0, 0.02, 0.04. Comparing with (|5.17|) we see 
that the quark mass corrections to the kinetic terms strongly influence the determination of 
Zh/Z m . We conclude that 1 — M^- can not be treated as a very small number. The difference 
between Z^ and Z m has to be included for any systematic discussion of the scalars within the 
linear sigma model! The neglect of this difference in earlier works [0- 14 1 partially explains 
the quantitative differences with our results. 

15 The exception from a systematic procedure of keeping only terms in the effective action for scalars that 
contribute to linear order in A is motivated by the well identified mechanism that induces a sizeable C,hV 2 
(cf . appendix |b|) . 

1 A similar observation holds in the pseudoscalar sector for m 2 g and we obtain m g 



393 MeV. 
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u h v 



Figure 5: The plot shows M ao as a function of ZJhV for Z p /Z m = 0.9, uj m v = uJ m v = —0.20 and three 
values of C,h.v 2 = (solid lines), 0.02 (dashed lines) and 0.04 (dotted lines). The upper curve in each band 
corresponds to CTq A3 = — A2/4 and the lower one to ctqXs = Aa/4. 



We are now ready to reexamine the mass splitting in the scalar octet ( 10. 2|) . Including 
the corrections arising from ( |11.3| ) we find 



3j 6 v - \m\u h v 



M l K% - Ml = 7_ 2 . (11.7) 

One sees that large negative values of TDhV can considerably increase the mass difference 
between the Kq and the ao for given ^v. The same holds for ChV 2 > 0. This weakens the 
argument of the preceding section against the association of the ao meson with the resonance 
ao(980). In fig. |5] we plot M ao as a function of To^v for Z p = Z m and three different values of 
ChV 2 , with bands corresponding to ranges of A 3 between <TqA 3 = — A 2 /4 (lower curves) and 
^0^3 = ^2/4 (upper curves). Here we have used ( |5.18| ) together with 



76 _ 1 Mh ~ K + (\M 2 Kt + Ml) 
72 2 Ml-Ml + uj m v (\M 2 K + Ml) 



(11.8) 



In fact, large quark mass corrections to the scalar kinetic terms seem the only plausible 
possibility for the choice of the ao(980). This would indicate a large mixing between two- 
kaon states and the ao(980) (see appendix |C|) . Consequently, it could explain why the cto(980) 
behaves in many respects similarly to a qqqq state even though it may belong to an octet of 
qq states. We should also point out that for large values of \uJhv\ the quark mass expansion 
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becomes questionable in the scalar sector. From the linearized expressions 



Mi 



m h- 2 (l6 V ~ \ m \^hV 



M 2 K% = ml+fav-~mlu h v > ) (11.9) 
we infer the ratio of the first order correction for M% as compared to m\ 



•y 6 v - \m\u) h v ) 9 Mi, - M' 2 

' - - ~ 0.10(0.43) , -0.03(0.33), -0.20(0.19) (11.10) 



mi 3 m( 

for QhV 2 = 0, 0.02, 0.04. Apparently, for (hV 2 = a good convergence of the quark mass 
expansion is only realized for the assignment ao(1320). For larger values of (hV 2 as infered 
from the leading mixing approximation in appendix [B] a reasonable convergence can also be 
obtained for a (980). 

Summarizing the various aspects of the problem of the correct assignment of the isotriplet 
belonging to the same octet as the _K"q(1430) we may state that the association a (1320) would 
make the understanding easier only in case of a standard kinetic term for h. Taking into 
account nonminimal kinetic terms there is no conclusive argument to rule out the ao(980) 
as a member of the scalar octet. For the latter assignment one expects important mixing 
effects with two-kaon states. Actually, such large mixing effects concern presumably only 
the ao and not the other members of the scalar octet. It may therefore be preferable not to 
include these mixing effects into the parameter UhV appearing in ( |11.3|) but to treat them as 



additional corrections to the ao propagator only. In this case the size of UhV is expected to 
remain small, \uh.v\ ^ 0.1, but the physical mass of the ao is related to M ao by an unknown 
factor reflecting the mixing. The value of M ao appearing in formulae like ( |11.5| ) or ( |11.11 ) 



below should then be replaced by an effective mass M ao somewhat above 1 GeV (say around 
1100 MeV). A natural mechanism of "threshold mass shifting" leading to a physical mass 
M ao = 980 MeV is described in appendix |C|. 



We finally want to show that the Gell-Mann-Okubo type mass relation ( 2.22|) is not 



2 

affected by linear quark mass corrections to the kinetic terms. Inserting M% o = M ao Z~^, 



M\* = M K *Z K l, M? = M f g Zfg into (2.22) one obtains to linear order in UhV the relation 



M 2 h = -A - V - \w h v (M 2 KS - M 2 ao ) . (11.11) 



The correction to the relation fl2.22|) is indeed only quadratic in A. For given M ao = 
1320 MeV and M K * = 1430 MeV the symmetry relation (|11.11|) yields for Q h v 2 = 

M /g ~ 1465 MeV. (11.12) 

On the other hand, for ( h v 2 = 0.02 and M ao = 1100 MeV we find 

M /8 ~ 1354 MeV. (11.13) 
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Figure 6: The plot shows Mf as a function of M ao according to the scalar Gell-Mann-Okubo relation 
Q2.22D with ( |11.3| ), Z p /Z m = 0.9, uj m v = uJ m v = -0.20 and fixed M K * = 1430 MeV. The bands correspond 
to values of A3CT0 between — A2/4 (upper curves) and A2/4 (lower curves) and we give results for (^v 2 = 
(solid lines), 0.02 (dashed lines) and 0.04 (dotted lines). 



Taking into account the uncertainties from the mixing with the scalar singlet s both values 
are consistent with the observed broad resonance^ /o(1300). In fig. |6] we have plotted 
Mf s as a functionP] of M ao in order to demonstrate the relative insensitivity of Mf g on the 
identification of the ao meson. For this purpose we have used (h = 0, 0.02 and 0.04 and 
\3<7q = ±\ 2 /4:. For each set of parameters we have determined UhV from fig. [|. Figure || 
demonstrates that for ChV 2 ^ 0.02 values of Mf 8 below 1500 MeV are preferred. On the other 
hand, a value of Mf$ below 1100 MeV would require a very substantial mixing of Kq with 
a state in d^py ((hV 2 ^ 0.04.). The identification of either the / (1590) or the / (980) with 
the octet seems therefore disfavored. 



In this context it is perhaps interesting to note that the branching ratio |l5j 



T (f (1300) -»• 2tt) , 
R= / J0K =v~12.5 (11.14) 

r (/o(i30o) - kk) 

is consistent with an octet assignment of the /o(1300). For a pure octet this ratio should be 
around 6 and a relatively small admixture of a singlet state could easily explain a further 
enhancement. (Whereas a pure singlet state would lead to a much smaller ratio R ~ 1.5 a 

17 The mass of the /o(1300) is not determined very precisely. It could easily be around 1400 MeV. 
18 In figs. H and ^| we have not distinguished between M 0o and M ao . Taking into account the additional 
mixing with two-kaon states for the ap the relevant axis actually shows M ao . 
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second solution for large R corresponds to a large octet-singlet mixing angle, i.e. tan$ s ~ 
—0.9 for equal sm 2 and hm 2 couplings.) On the other hand, the branching ratios of the 
higher mass resonance /o(1590) seem compatible with a singlet with large cubic coupling 
~ p 2 s, but not with an octet. For the /o(980) an assignment is difficult in view of the 
presumably large mixing with two-kaon states. 

In summary, two natural scenarios for the scalar nonet seem to be compatible with 
the parameters of the linear meson model extracted from the pseudoscalar sector: For one 
scenario the isotriplet corresponds to the ao(980) and the singlet (or dominant ss state in 
case of large \$ s \) is associated with the / (980). These four states are largely influenced by 
mixing with two-kaon states. The other f state of the nonet (dominantly an octet state 
in case of small |$ s |) corresponds to the /o(1300). In this case a relatively large mixing in 
the strange sector with d^py (large C,h,v 2 ) explains why the Kq(1430) has the highest mass 
in the nonet. The other scenario has a larger average mass of the octet. The triplet is 
associated with the proposed a (1320) which is not far below the doublets _K"g(1430). The 
ao(980) and /o(980) are dominantly four-quark states or KK molecules in this case. Again, 
the octet state is the /o(1300). The singlet corresponds^ either to /o(1590) or its width is 
too large to be detected. Reliable information about the value of Ai would certainly be of 
great interest for further pinning down the possible options. 



12 Higher derivative contributions 

In this section we investigate deviations of the meson propagator from the approximated 
form Gi = (Ziq 2 + M^) -1 . We will first concentrate on the flavored pseudoscalars. In the 
language of sect. |6| we want to make an estimate of the corrections AZi ( |6.13| ). Within a 
systematic quark mass expansion we need to order A 2 only the q 4 correction to the inverse 
propagator in an approximation where it is independent of the quark masses. This correction 
arises from a term involving four derivatives 

£ kin (4) = -j 1 Tr (d 2 md 2 m - qld^md^m) (12.1) 

where q 2 is chosen according to ( |6.19| ). It involves one additional parameter H m which 
determines the ratio 

u m P 2— M 2 -M 2 , 
— = 1 + ^ = 1- -H m — * * . (12.2) 

As a first observation we notice that H m receives contributions from the mixing of the 
pseudoscalar mesons with other states. We infer from appendix that the mixing with 
other 0~ + octets indeed induces higher derivative corrections because of the momentum 
dependence of the propagators of the additional states that are integrated out. If we assume 
that the dominant contribution to uo m arises through mixing with other states we can identify 

19 There are other not so well established resonances /o(1510) and /o(1525) Jh| which may be identical 
with the /o(1590) or also be possible candidates for the singlet. 
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in ( |C.7|) oj m with u$ and find 

where we used M P = 2280 MeV(2670 MeV) (cf. appendix and table | in appendix g). 
One may also estimate K 8 according to ( |6.19D and finds 

K { s p) ~ 0.002(0.001) . (12.4) 

This is indeed negligible for the wave function renormalizations as compared to uJ m v- 

A different contribution to the higher derivative term (|12.1|) arises from loops of meson 



fluctuations. For an estimate of their importance we use the modified loop expansion of fi~6 | 
("systematically resummed perturbation theory"). This allows to compute the deviations 
of the inverse propagator from the form q 2 + M 2 in terms of Mf and effective 1PI cubic 
vertices. It is crucial in our context that instead of "classical vertices" only the 1PI Green 
functions appear in the perturbative series since only the latter are directly calculable from 
the present phenomenological analysis. Also the loop expansion is only used for "higher 
order couplings" where it converges reasonably well. (We do not expect a good convergence 
for quantities like Z h /Z m etc.) Let us write the inverse propagator for a member of the 
pseudoscalar octet as 

G-\q 2 ) = q 2 + M 2 + tM 2 ). (12.5) 

Here we have subtracted from the usual self energy Sj(g 2 ) those parts which are already 
contained in the effective wave function renormalizations Zj and mass terms M i 

Uq 2 ) = £,(g 2 ) - £,(0) - £ fefoS) - Si(0)) . (12.6) 
This definition implies 

Ei(0) = t^ql) = (12.7) 

and we use g 2 = — rr? m . In consequence, Sj(g 2 ) contains only contributions to higher deriva- 
tive terms. The dominant one-loop contribution to S m (g 2 ) for the pseudoscalar octet is 
depicted in fig. [7|. It involves the propagation of a scalar and a pseudoscalar in the loop 
and we therefore need the cubic coupling 72 (cf. section ||). We observe that the subtrac- 



tion ( 12.6 ) makes the usual one-loop expression ultraviolet finite. In fact, if the momentum 
dependence of the effective three-point vertex is not too strong the momentum integral for 
the difference S m (g 2 ) — S m (0) is dominated by momenta in the range between the masses of 
the two particles propagating in the loop. We will use here the approximation of a constant 
cubic vertex which we approximate by its value to lowest order in the quark mass expansion 
given by 72. We also neglect the octet mass splitting for the particles propagating in the 
loop for which we use average masses m m = 412 MeV and rrih = 1394 MeV. To lowest order 
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mi(g) 



Figure 7: One-loop diagram for the dominant contribution to Ei(g 2 ). 



in the quark mass expansion we are interested in the effective coupling H m (|6.12 ) which is 
given by 



Hr 



_B_ (ZrA 

d( i 2 v ^ / 9 2 =? 2 



;i2.s) 



We are interested in the momentum range — q 2 < (m^ — m m ) 2 for which the one-loop 
contribution is given by 



^V)-^(o) 



48tt 2 
1 

Q 

x In 



m m , m m + m h 



In 



12.9) 



m m + m h ) 2 + g 2 + J (m m - m h ) 2 + q 2 ' 



^{m m + m h ) 2 + q 2 - ^ (m m - m h ) 2 + q 2 J 

Here we have neglected contributions ~ 7i, t| as well as the 77-77' mixing. For q 2 = —m 
this yields the one-loop contribution to H m 



H 



(i) 



m 2 m (m 2 h -m 2 m ) m m 



48vr 2 m 4 



; 12.10) 



+ m hX / m 2 h - 4m 2 m — + 



y/m h + 2m m + ^m h - 2m r 
m li ' m \~ ^ m mJ ~* + 2m m - - 2m r 



In 



Using 72 ~ 8000 MeV and a; m ?; ~ —0.20 (see sect. |13|) this results in 



if^ ~ 1.64 ■ 10 MeV" 



'12.11^ 



or 



~ 0.012 . 



12.12) 
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We see that loop corrections to the higher derivative terms are negligible as compared to 
contributions arising through the mixing with other states. We will therefore assume that 
the higher derivative terms are dominated by such mixings and estimate 

^ = ^^-0.95. (12.13) 

Urn l + 5 ( f 

Even though (p m — u m )/uj m is formally not suppressed by powers of A we see that this 
higher derivative effect is actually small. 

In the 77-77' sector we need information about the momentum dependence of z 8 (q 2 ), z p (q 2 ) 
and co(q 2 ) (|6.26|) . The relevant quantities are d 8 = z 8 (—M 2 ,) — l,d p — z p (—M 2 ,) — 1 and d w = 
(£>(— M3) — l))/Cj where we remind the reader of the definitions z$(— M 2 ) = z p (—M 2 ) = 1 
and u(—M 2 ) = Cj. Since the mass difference M 2 , — M 2 exceeds substantially M\ — M 2 the 
higher derivative corrections in the 77-77' sector could be somewhat larger than for the flavored 
mesons. Altogether, the contributions beyond the lowest order in the derivative expansion 
involve four additional dimensionless parameters, 5^, d 8 , d p and d^. Their absolute size 
is expected to be small if the derivative expansion converges. Since the predictions of the 
lowest order in the derivative expansion come already very close to the experimental values 
of f v and f v ' it seems not difficult to achieve agreement with observation by using small but 
otherwise arbitrary values for these four parameters. The expansion in powers of o"o may 
lead to some approximate relations between 5^, d 8 , d p and d u but we will not pursue this 
issue here further. 

Instead, we conclude this section by a description of the predictions of the "leading 
mixing approximation" . For this purpose we assume that all corrections to the kinetic terms 
- both, quark mass and higher derivative corrections — are due to a mixing with higher 
states contained in the divergence of the axialvector field d^p\. The formalism is described 
in appendices [B| and |C| All parameters appearing in the kinetic terms can be computed in 
terms of masses and couplings of the vector- and axialvector fields. Most of these couplings 
can be determined from observation (cf. appendix 0). There remains essentially only one 
important free parameter Z P which appears in the term (1/4) Tr ^Zp(d fJi py) 2 + Zp{d lx p^ 2 \. 
This parameter determines the strength of the higher derivative terms induced by the mixing 
with M 2 P ~ Zp 1 in In fig. |a we plot the values of M v /M° xp , f v /f° xp and frf/fif* 

as functions of Zp. (The higher derivative corrections vanish for Zp = 0.) For this plot we 
use x p = 1 (cf. appendix H) and employ the leading mixing approximation which assumes 
1$ and u m /u m given by ( 12.3 ), (|12.12| ). The higher derivative contributions in 



( |6.26|) and (|A.26| ) are now included according to ( |C.9| ), ( |C.13| ) and (|C.14j ). We see that a 



reasonable picture emerges for Zp ~ 0.22. It is consistent with the assumption that the 
nonminimal kinetic terms for the pseudoscalars are dominated by the mixing with d^p^X or 
the "partial Higgs effect" . The remaining differences of the curves from one can reasonably 
be attributed to subleading effects beyond the leading mixing approximation, as described 
in the more general framework of the main text. For a demonstration we also show in 
fig. [8]b the situation which arises if in addition to the contributions from the leading mixing 
approximation one also includes a nonvanishing U v in the kinetic terms (f4.6|) . For this plot 
we have chosen U 9 such that Z p /Z m = 0.9 (cf. flS.10|) . The agreement with observation 
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Figure 8: The plots show the curves for M v /M^ (solid line), f v /f° xp (dotted line) and f n >/f°? p (dashed 
line) in the "leading mixing approximation" as functions of Zp for x p = 1, Z p /Z m = 1 (a) and Z p /Z m = 0.9 
(b). 
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M K ± 


M K o 


M v > 


u 




135.1 


492.4 


497.7 


957.8 


92.4 


113.0 



Table 1: This table shows the phenomenological input used in this work. All values are given in MeV. 
The charged meson masses are electromagnetically corrected. 



improves and the optimal value of Zp is shifted to somewhat smaller values. We emphasize 
that the leading mixing approximation is complementary to the formal expansion in powers 
of quark masses. It is encouraging that a simple mechanism (the partial Higgs effect) can 
apparently explain the dominant parts of the parameters appearing in the systematic quark 
mass expansion. 



13 Results 

For the convenience of the reader we summarize in this section the results of a numerical 
solution of our equations. We observe that the pseudoscalar sector can be treated indepen- 
dently from the scalar sector. For the flavored pseudoscalars there are two small parameters 
whose influence is rather modest, namely Z p /Z m — 1 and 5 W = ZJ m /uj m — 1. Three additional 
small parameters d w , dg, d p characterize the most general form of the higher derivative terms 
in the rj-rf sector. For the first two lines (a) and (b) in our tables we use the first order in the 
derivative expansion, i.e. 5^ = d^ = d$ = d p = 0. We present two values Z p /Z m = 1.0 and 
0.9. The value of the dominant free parameter u m v is chosen such that M v comes out close 
to its experimental value (cf. fig. [I]). Going beyond the first order in the derivative expan- 
sion we include in line (c) the higher derivative corrections corresponding to a nonvanishing 
5 U . It is taken according to the leading mixing estimate ( |12.13j ) such that u> m /uJ m = 0.9. 



Finally, we present the results of the leading mixing approximation in line (d). Here all 
quark mass corrections to the kinetic and all higher derivative terms are determined from 
the simple assumption that they are induced by the exchange of the axialvector field d^p^ A . 
Line (d) corresponds to fig. [8|b with Zp ~ 0.16. In the leading mixing approximation u m v 
and uJm/uJm are not anymore free parameters and can therefore not be adapted to fix M v to 
its experimental value. For line (d) the values of M v , f v and differ somewhat from those 
obtained for the optimal value u m v = —0.20. The leading mixing approximation comes 
nevertheless quite close to the experimental results. In table |l| we give our input values for 
the pseudoscalar sector. Table |2| shows the four different combinations of the parameters 
u> m v, Z p /Z m and uj m /uj m used for our numerical analysis. The numbers given in the text of 
this paper correspond to line (b) which may be considered as our best values to first order 
in the derivative expansion. 

The first step is to solve for the Z-factors and to extract the values of the zero 

momentum parameters M i , f { . The corresponding relations are ([2.6|), (|5.6|) , (|6.3D, ( |6.10| ) and 



the results are found in tables |3| and |j. We should point out the sizeable differences between 



r2 



M/ and M i and similarly for the decay constants. They are due to the large value of \oj w v\ 
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UJ m V 






ZJ m v 


(a) 


-0.22 


1.0 


1.0 


-0.22 


(b) 


-0.20 


0.9 


1.0 


-0.20 


(c) 


-0.22 


0.9 


0.9 


-0.24 


(d) 


-0.17 


0.9 


0.95 


-0.18 



Table 2: The table shows the four different combinations of the parameters uj m v, Z p /Z m and ui m /ZJ m 
used in sect. [13[ Line (d) corresponds to the leading mixing approximation for which uj m v and uj m /uj m are 
computed and do therefore not play the role of free input parameters. 





Ztt 


Z K ± 


Z K ± — Z K o 


M ± 
MeV 


M K± 

MeV 


M k0 
MeV 


(a) 


1.22 


0.89 


0.0063 


149.2 


465.3 


468.7 


(b) 


1.20 


0.90 


0.0058 


148.0 


467.9 


471.4 


(c) 


1.24 


0.88 


0.0068 


150.7 


462.2 


465.4 


(d) 


1.18 


0.91 


0.0052 


146.5 


470.9 


474.6 



Table 3: Values for wave function renormalizations Zi and zero momentum mass parameters Mi. 





U 

MeV 


/jf± 

MeV 


O"0 
MeV 


( Z m \ 1 v 
\ Z h ) MeV 


fv \1/ 2 
/ Z m \ 1 w 

V Z h ) McV 


(a) 


83.7 


119.6 


53.9 


24.2 


-0.69 


(b) 


84.3 


118.9 


53.8 


23.3 


-0.67 


(c) 


82.8 


120.4 


54.1 


25.3 


-0.71 


(d) 


85.2 


118.2 


53.7 


22.2 


-0.66 



Table 4: Expectation values of scalar singlet and octets. 



56 





m 9 

McV 2 


V 

MeV 


A 2 


UJ 


(z h y/ 2 72 

V Z m J MeV 


72" 
MeV 2 


72 w 
McV 2 


73" 

McV 2 


(a) 


(390.9) 2 


6814 


17.0 


0.35 


4780 


(340.1) 2 


-(57.4) 2 


-(80.9) 2 


(b) 


(392.9) 2 


6447 


21.3 


0.30 


4942 


(339.3) 2 


-(57.7) 2 


(46.1) 2 


(c) 


(388.6) 2 


5964 


17.5 


0.33 


4405 


(333.8) 2 


-(55.7) 2 


-(38.6) 2 


(d) 


(395.2) 2 


6011 


26.8 


0.25 


5167 


(338.8) 2 


-(58.2) 2 


(112.3) 2 



Table 5: Parameters of the linear meson model. 





M n 
McV 


fv 
MeV 


w 

MeV 


h 
U 


h 


h 

Jri 


V 


(a) 


550.8 


106.6 


93.5 


1.15 


1.01 


1.09 


1.25 


(b) 


546.9 


113.4 


83.8 


1.23 


0.91 


1.16 


1.12 


(c) 


549.1 


103.5 


88.8 


1.12 


0.96 


1.06 


1.19 


(d) 


536.6 


111.3 


82.9 


1.20 


0.90 


1.14 


1.11 



Table 6: "Predictions" for M v , f v and f n >. 



In particular, the difference f K — f n is almost twice the value of fx — fJ In the next step 
we determine in table [5] the parameters of the linear meson model from (|5.2j )- (|5.4j ), (|5.6p . 
( |8.12|) and ( |2.10| ), Q4.29Q . We observe that the infered value of the cubic coupling v depends 
only moderately on the details of the effective meson model. In contrast, the uncertainty 
for the quartic coupling X 2 remains substantial. Table [5] also contains information about the 
cubic coupling 72 between two pseudoscalar octets and the scalar octet as well as for the 
coupling 73 between the pseudoscalar octet, the pseudoscalar singlet and the scalar octet. 
Even though the sign of 73 remains undetermined we find 1 73/72 1 ^ 0.1. The coupling 72 
therefore dominates the decay of the ++ mesons. 

We are now ready to compute the mass matrix elements of the 77-77' sector using ( |5.4j ). 
The eigenvalues M v , M v > follow from (|6.29|) the mixing angles from (|6.27| ) and the relations 
for f n and f v t are given by ( A.23j ) and ( |9.9| ). These are our main "predictions" for observable 
quantities. They are displayed in table |6|. We find a very satisfactory agreement with 
experiment for line (b). The uncertainty in the "prediction" for f v and f v > is reflected in the 
differences as compared to lines (a) and (c). Also the leading mixing approximation (d) is 
not too far from experiment, even though contributions beyond this approximation need to 
be included. The general tendency of the higher derivative contributions in the 77-77' sector 
is an enhancement of f v / and a decrease in M v and f v (cf. fig. |8|). From table [7] we note that 
the octet decay constant f v s is rather close to fx whereas the corresponding singlet decay 
constant almost equals f n . In table [5] we also show the mixing angles in the 77-77' sector. 
We find a large mixing between 77 and 77' with an important dependence on the momentum. 
The mixing for q 2 = — M 2 , is substantially larger than that for q 2 = —M 2 . The mixing 
corresponding to p (rj) ~ —13.7° is somewhat smaller than earlier estimates from chiral 



57 





fvs 
MeV 


/t,'8 

MeV 


McV 


fri'O 

MeV 


o P (v) 




m v 
MeV 


~W, 

v' 


MeV 


(a) 


116.6 


132.1 


95.3 


107.9 


-15.6 


-31.2 


839.1 


0.77 


0.25 


(b) 


117.1 


124.2 


96.3 


102.1 


-13.7 


-28.0 


865.9 


0.82 


0.28 


(c) 


116.0 


126.6 


94.0 


102.6 


-16.3 


-31.1 


839.8 


0.77 


0.22 


(d) 


117.7 


120.4 


97.5 


99.8 


-14.7 


-28.6 


843.5 


0.78 


0.30 



Table 7: Mixing angles in the 77-77' system and isospin violation in the decay constants. 



LU h V 


ChV 2 


^0 


M K * 
McV 


Zh_ 

Z m 


m h 
MeV 


76 v 
MeV 2 


72 

MeV 


MeV 


M h 

MeV 


0.0 


0.00 


0.0 


1430 


0.35 


1407 


(254.7) 2 


8394 


1360 


1453 


-0.2 


0.00 


0.0 


1430 


0.39 


1335 


(241.6) 2 


7963 


1178 


1541 


0.0 


0.00 


9.4 


1430 


0.35 


1378 


(382.6) 2 


7973 


1267 


1480 


-0.2 


0.00 


9.4 


1430 


0.39 


1307 


(363.0) 2 


7564 


1097 


1570 


0.0 


0.02 


0.0 


1295 


0.42 


1274 


(230.6) 2 


7601 


1232 


1315 


0.0 


0.04 


0.0 


1144 


0.54 


1126 


(203.7) 2 


6715 


1088 


1162 



Table 8: Masses and parameters of the scalar octet. 



perturbation theory || where the size of the mixing angle was extracted indirectly from the 
requirement of a realistic value for M v . On the other hand, the mixing corresponding to 
9 p (rf) is larger. We should point out that our direct method of computing all elements of 
the matrix for the inverse propagator in the 77-77' system is quite different from the indirect 
consistency requirement for M v . We furthermore see in table a large deviation of m 2 /M 2 , 
from one despite the fact that this difference is formally only a quadratic term in the quark 
mass expansion. We also present in table [7] the isospin violation in the decay constants 
(16.241) , f K o — f K ±. It is reduced significantly as compared to the value obtained to lowest 



order in the quark mass expansion Q2.12j ) or in chiral perturbation theory [0] . 

Our results in the scalar sector depend in addition on Z h /Z m , A 3 cr 2 , ZJ h v and (h^ 2 - We 
use A3CT 2 , UhV and (^v 2 as input parameters together with a fixed value M K * = 1430 MeV. 
From (|4.33| ), (|4.34j) we determine m\ and j^v as functions of these three couplings. Here we 



use our "optimal values" for A 2 , u, etc. corresponding to the second line (b) in tables [2|-|7|. 
We determine Z h /Z m as a function of A3, oJhV and (hV 2 according to 



Zh_ = m 2 + ct (3ct A 2 + v) + |(9ct A 2 -v + 12ct q 3 A 3 )[/^ - f vl , | 



J in. 



1 + \uj h v 



M 2 K , 



The results of this analysis and, in particular, the masses of the lowest lying ++ octet are 
given for several values of A3, W^v and C,hV 2 in table ||. Here A 3 cr 2 = 9.4 corresponds to the 
"maximal value" A 2 /2 compatible with a convergent expansion in o"o (cf. sect. |10f) . From an 
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estimate of (hV 2 due to the exchange of the vector field d^py (see appendix |B|) we learn that 
the two last lines in table |8| are preferred. Additional large mixing effects (see section |ll| and 
appendix 0) may further lower M ao and lead to a mass of the isotriplet ao consistent with 
the observed resonance a (980). The scalar partner of the 77 appears to be associated with 
the broad resonance /o(1300). We finally emphasize the large deviation of Zh/Z m from one 
which underlines the importance of nonminimal kinetic terms in the linear meson model. 



14 Conclusions 

In this paper we have investigated the effective action for the linear meson model. Including 
the discussion of the vector and axialvector fields from appendix |B| a fairly simple picture 
emerges. Expressed in terms of scalar fields $ and vector fields py, p\ the quark mass 
expansion seems to converge rather well for the three light flavors. The same holds for 
the derivative expansion, leading to an approximate momentum dependence of propagators 
~ [Zq 2 + M ) . The only exception from this picture seems to be the scalar isotriplet 
ao (980) which can be explained (see appendix 0) by a large contribution of two-kaon or 
four-quark states. 

The divergence d^p'X has the same quantum numbers as the pseudoscalar octet plus 
singlet. We have estimated the resulting mixing effect or, equivalently, the terms induced 
in the effective action for $ from integrating out d^p^ ("partial Higgs effect"). We find 
a large non-minimal kinetic term which induces substantial quark mass corrections to the 
kinetic terms for the pseudoscalars (u)!fi)v ~ —0.15). This effect remains compatible with a 
converging quark mass expansion for the flavored pseudoscalars 7r and K. On the other hand, 
an investigation of the masses M v and as well as the decay constants f v and f v > shows 
that the non-minimal kinetic term induces in turn a large momentum dependent mixing 
in the 77-77' sector. Here we find that contributions which are of third or higher order in a 
formal quark mass expansion are comparable in size to the contributions arising to second 
order. It is the nonlinearity generated by this large mixing which leads to "predictions" for 
M v , f v and f v ' which are compatible with experimental observations. This explains why the 
measured value of f v is quite far away from its value for zero quark masses^. It is amazing 
to see that the nonlinearities in M v , f v and f v > as functions of the variable u m v all conspire 
such that a common value of u m v can explain simultaneously these three quantities. Even 
more, our phenomenological estimate u m v ~ —0.2 is quite close to the estimate from the 
partial Higgs effect, ~ —0.15. The latter involves completely different quantities like 

the pirn coupling and the masses of the axialvectors! 

It will be very interesting to see if our "phenomenological" picture of the effective linear 
meson model can be obtained from the solution of a flow equation similar to ||. It would be 
highly nontrivial if the couplings would come out in such an approach in a range consistent 
with the analysis of the present paper. This concerns, in particular, the quartic coupling 

20 We emphasize once again that the nonlinearities in A or the poor convergence of an expansion in A 
appear only in the eigenvalues M v and M v >, the mixing angles 6 p {rj), P {r]') and the decay constants f v , f v >. 
The matrix elements of the 77—77' propagator converge satisfactorily. 
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A2 which was found in || to be essentially determined by an infrared fixed point behavior. 
Furthermore, a more systematic analysis of the (axial)vector meson sector including quark 
mass effects should lead to a quantitative estimate of several effective cubic and quartic 
vertices relevant for the decay properties of these mesons. Beyond a successful explanation 
of the observed values for f v and f n > our results constitute the "phenomenological basis" 
for interesting further developments. They also shed light on the important question of the 
convergence of the quark mass expansion. 



Appendices 

A Meson decay constants 

In this appendix we discuss the meson decay constants within the linear sigma model. Most 
results displayed here are well known from current algebra and are simply rephrased in a 
somewhat different language. The only slightly delicate issue concerns the choice of the 
normalization of fields. This determines the appropriate definition of wave function renor- 
malization constants. A careful treatment of these constants is relevant for quantitative 
relations between decay constants and meson masses as discussed in the main text. We 
adopt here definitions of f n , fx, f v , etc. which are directly related to measured partial decay 
widths of the corresponding mesons. 

The weak leptonic decay of the charged pion, ix~ — > involves the effective three 

point vertex (7 being the Euclidean analog of 75) 



r - 



d 4 p n d 4 p fl d^pu 
(2tt) 4 (2tt) 4 (2tt) 4 



ff 2 cos?9 e 



™ i (A.l) 

x 7r~(p 7r )/Z(p M )7p(l + 7)z/^(-p F )(27r) 4 5(p 7r - p^ - p?) + h.c. 

Here we have projected onto the leptonic V — A structure following from virtual W -exchange 
with M\y the py-boson mass, g the weak gauge coupling and $ c the Cabibbo angle. (Analo- 
gously, the effective vertex for the charged kaon decay is obtained from (|A.1|) by the replace- 
ments F£ — * Fx, 7r~ — > K~ and cos$ c — ► sini? c .) The vertex function F% can depend only 
on two independent momenta (e.g., p n and pv) and the leptonic pion decay involves its value 
for on-shell momenta, p\ = —M%±, p 2 ^ = —M 2 , p§ = 0. In the present context we neglect 
the dependence of F£ on the leptonic momenta and use the parameterization^] 

F! = pPFM)- (A.2) 
Our task is therefore the evaluation of the pion decay constant 

f* = FM = ~Ml ± ) (A.3) 

21 A term ~ p~ would not contribute to the pion decay anyhow, since on shell we have i>vv{— , p^) = 0. 
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which is determined experimentally from the leptonic decay width of the pion (up to elec- 
tromagnetic corrections) 

iw = ^K^-fl (i - wij^ 2 ^ ( A - 4 ) 

derived from flAl|) with G F = \/2# 2 /(8M^). 

In order to compute F w we consider the linear a-model coupled to external currents. For 
this purpose we replace all derivatives acting on $ by covariant ones 

D»$ = « - % -\ z R^ + ~$A*2£ = - % -V? [X z , d>] - % -A* {A„ $} . (A.5) 

Here the vector- and axialvector currents and are related to the left and right handed 
currents and R%, respectively, by 

Ut = V?-A>i 

(A.6) 

Rt = V? + A'i . 



By this replacement Tkm = / d 4 x£kin (|4.6| ) becomes a functional of the (local) background 
fields V 11 and A* 1 . Current conservation is automatically embodied in this construction. The 
coupling of mesons to V^-bosons can be infered via the identification 

1,2 ' (A.7) 

L£ 5 = gsm# c W£ 2 . 

Once the couplings of mesons to VF-bosons are known the couplings to lepton pairs follow 
by insertion of the field equation 

W<P ) = -f ^Sjp / - ^1- (1 + WW • (A.8) 

Here ip stands for the lepton doublets and we will neglect p 2 as compared to Mjy. Extracting 
the coefficient linear in L^ 2 



(2vr) 4 

d 4 p 



(A.9) 



2V2J (2tt 

we can relate F n to the part of K£ 1 which is linear in 7i"i 



7r-(p)F£(p) [L^{p] - zL 2 »] + . . . 



Kl 1 (p)= % -p^F 7T (p)n 1 (p) + ... . (A.10) 
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The discussion of leptonic decays of charged kaons is analogous with 



KM 



-yF K {p)K A { P ) + 



(A.ll) 



Here 711,2 and are the pion and kaon fields, 7r = -^(vri + ^2), K = -^{K± + iK 5 ), 
with standard normalization of their kinetic terms such that their inverse propagator (two 
point function) in the vicinity of its zero at timelike momenta is given by p 2 + M^±(Mjl±). 

For an evaluation of the terms linear in the currents V^, At, and linear in the meson 
fields one covariant derivative should be replaced by 



D»* — -\v?[\ Mi (<&>] - \A»,{\ Mi ($))} 



(A.12) 



with ($) = diag(<y9 u , ifd, <p s ) the expectation value of $. All other covariant derivatives have 
to act as simple derivatives on the meson fields. Consider first the limit of vanishing quark 
masses where ip u = tpa = ip s = Wq. The term linear in the vector current vanishes whereas 
the axial current couples linearly to the pion field as a consequence of spontaneous chiral 
symmetry breaking, i.e. -D M $ — > — ia A^X z . The relevant term in £ kin fl4.6|) reads 



kin 



and, with Z n 



Z^ + X (p a 2 ) + Upao (cf. sect. 



U v a ) Z-'aoA^Tr {X^d^m} 
H), one infers 



(A.13) 



--<j Tr {A 2 ,<9 M m} 



(A.14) 



In the limit of vanishing quark masses the fields m z are already properly normalized and we 
can identify tti = mi, K\ = 777.4. This yields 



f n = 2cr 



(A.15) 



The proportionality between f n and <Jq is no surprise. It is well known that in the limit of 
vanishing quark masses the pion decay is related to the non-conservation of the axial current 
induced by chiral symmetry breaking. 

Going beyond the limit of vanishing quark masses the expectation values of (p u , (fa and 
ip s are different. Nevertheless, there is no term linear both in the vector current V£ and 
a pseudoscalar meson. This follows from the discrete symmetries C, V and is related to 
the observation that the first term on the right hand side of ( [A. 12 ) is hermitean whereas 
the second is anti-hermitean. On the other hand, the vector current couples linearly to the 
scalars as, for example (see eqs. ( fl4|) , (|4.11|) ) 



Tr (D^y 



--AtTt 
2 z 



-^K M Tr 

2 z 



[A„ (*>] (z h 1 



%h + -j=Z s 



2 d fl s 



(A.16) 
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For our computation of f w we can omit the term linear in Vj*. 

Let us consider first the approximation where the kinetic term ([Of ) is truncated to 
Z v Tr (D^Q)^ (D^). In this limit m describes already the properly normalized pion and 
kaon fields. One finds 

Ki = \f*i(p)Z v Z^ { Vu + <p d ) (A. 17) 

or 

U = Z v Z m '(<p» + <p d ) ■ (A. 18) 

Similarly, we find for the leptonic decay of the charged kaons 

fK± = Z ip Zn*(<p u + <p a ) (A. 19) 

and define 

f K0 = Z V Z^ (<p d + <p a ) . (A.20) 

Noting Wq — | (<p u + ifd + ip s ) and observing that in this approximation Z m = Z v one arrives 
at Q 

fn + fi<± + f K o = 6a z| = 6a . (A.21) 

With the definitions A u = Zl^u ~ etc. and (h) = ^piz\ ($ a ) = 2z\ (($) -W ) = 
2 (Zh/Zm) 1 diag(A u , A d , A s ) we obtain the relations (|2.6|). 

Next we consider the more general kinetic term (|4.6|) . The first effect is a nontrivial wave 
function renormalization between the fields mi and 7Ti, i.e. 

_ i _ i 

mi = Zir 2 ni , m 4 = Z K \K^ . (A. 22) 

_ i 

This effect multiplies by a factor Ztt 2 and similarly for f K ±, f K o. Here we note that Z n , 
Zx± and Zj^o should be defined at the corresponding poles such that the inverse renormalized 
two-point function is approximated in the vicinity of the pole by q 2 + M 2 ± with M n ± ~ 
135.1 MeV the physical pion mass (after subtraction of electromagnetic effects). The second 
effect reflects the modification of the general kinetic term into which ( |A.5| ) is inserted. Since 
the axialvector field A^(q) is needed for on-shell momenta, we conclude that all kinetic 
terms must be evaluated at the poles. Expanding the inverse propagators around q 2 = 
—M 2 , knowledge of the coefficient linear in q 2 suffices for a computation of the meson decay 
constants. For the pions this is given by Z v and the evaluation of the full kinetic term gives 
a factor Z n in the formula for /„.. This can be seen directly by inserting ( A.5|) into the 



contributions ~ X~ U v (|8.6| ), (|8.9| ) and using the relations (|8.10| ). In summary, the total 
effect of the generalized kinetic term is a multiplication of f w with Z\l 2 . Similarly, fx± and 
fxo are proportional to Z^l and Z^o, respectively. This explains the relations (|6.9|) . We 
emphasize that Z n , Z^± and Z K o should be evaluated from the coefficient linear in q 2 in an 
expansion of the inverse propagator around q 2 = —M 2 ± , — M|- ± , — M|- , respectively. More 
precisely, they are defined by ( |6.11| ) with q 2 . replaced by —M 2 . 



We finally extend the discussion to the decay constants of the non-flavored pseudoscalars 
f n o, f v and f v >. By eqs. (|9~^ ) and (|9.10|) they are related to the couplings of these mesons 
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to the corresponding components of the axialvector currents, or in a different language, the 
expectation value of the corresponding current between the vacuum and the meson state. 
For instance, f v8 parameterizes the coupling of rj to the current A 8 whereas f v > corresponds 
to the coupling of the rf to the singlet current. For a comparison of f v8 with f n we therefore 
have to replace Tr ({Ai, (<3>)}Ai) Z}J 2 by Tr ({As, ($)}Ag) Z\j 2 . The ratios of /„-o, f v8 and f v >o 
to /„- are then easily computed 



i 



f v8 / Z 8 \ 2 p u + ip d + 4ip s f Z 8 \ 2 2f K ± + 2f K0 - f n 



fn \ZirJ 3(v?« + y?d) \Z W J 3f w 

f v o _ {Z 8 y 2(tp u + tp d + tp s ) _ /z 8 y J K ± +J K o + / w 



(A.23) 
(A.24) 
(A.25) 



Similarly, the octet and singlet decay constants for the normalization of the r]', f v ' 8 and f^Q, 
are given by 

^ = ^=(^YU-M^. (A.26) 

Here z p ( — M 2 ,) accounts for higher derivative effects and is normalized for q 2 = —M 2 accord- 
ing to z p (—M 2 ) = 1. If the higher derivative effects are omitted one also has z p {— M 2 ,) = 1. 
Furthermore, if we neglect the mixing effects (or for z 8 (q 2 ) = z p (q 2 )) we can identify z p (q 2 ) 
with z p (q 2 ) appearing in ( |6.26| ). We note that ( [A. 26 ) is appropriate for a definition of Z p at 



q 2 = —M 2 . If one instead would define Z p at q 2 = —M 2 , the factor z p (— M 2 ,) 1 ^ 2 would be 
absorbed by this alternative definition. 



B Vector mesons 

In this section we briefly discuss^ the vector and pseudovector fields and their interactions 
with scalars and pseudoscalars. This will permit us to estimate the part of the effective 
interactions involving four (pseudo) scalars which is induced by the exchange of vector fields. 
We introduce the fields and p R as hermitian 3x3 matrices which transform as 8 © 1 
under SUi(3) and SUr(3), respectively, being neutral with respect to the other part of the 
flavor group. The interaction with (constituent) quarks 

£pqq = ^ Z q9pqq {IlI^PW + QrI^PrCIr) (B.l) 

22 See also and references therein. 
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respects SUl(3) x SUr(3) and is also consistent with left-right symmetry (pi <-> Pr) and 
charge conjugation (p L — >■ — pj, p# — >■ —pE)- The invariant kinetic and mass terms read 

£ p , 2 = ^Tr(^p L ,-^p iM )(^-9X) + ^Tr(^) 2 

Z 1 i *" (R2) 

+ ^ Tr p£) 2 + -roj Tr p L „p£ + -m p 2 Tr p£ Tr p iM + (L - i?) 

where ^ 

Pl,r = Pl,r ~ 3 Tr Pl,/? = ( B - 3 ) 

denotes the octets and Tr p^ L R represents the singlets. In general, the field p M can de- 
scribe spin-one and spin-zero (~ d^p^) particles. For a p , a' — > the spin-zero components 
decouple and the fields p M only describes spin-one bosons, = 0). In the opposite limit 

a p , a' p — > oo there remains a kinetic term only for the spin-one bosons whereas d p p^ can be 
determined algebraically from the field equations. 

There is only one possible cubic coupling involving two scalars or pseudoscalars and the 
vector octet. To lowest order in a derivative expansion it reads 

= A^g^J Tr [(3„$t$ _ $t^$) f L + (q^\ _ $^$t) pM] (B 4) 

The appropriate value of Z will be determined later such that g pn7T appears in the width of 
the decay p — ► 7r7r according to 



9 2 (m 2 - 4M, 



3 

2\ 2 



Using the experimental values T(p — > 7r7r) ~ 150 MeV and M p ~ 770 MeV this yields 

Q p7T7T ^ 6.0 . (B.6) 

A second coupling appears for the singlets 

C*> p = ^ZTr fa&Q - (Tr - Tr p£) . (B.7) 

We note that similarly the couplings of the quarks to the vector and pseudovector singlets 
could be different from the octet couplings leading to a modification of ( B.l|) . In the following 
we will neglect for simplicity the differences between the singlets and octets (i.e., g' p7Tn = 
g P Tm, rn'p = m 2 p , a' p = a p ). Accordingly, we give for the quartic couplings involving two 
(pseudo) scalars and two (pseudo) vectors only those appearing for the octets and extend 
them to the singlets. In contrast to the cubic coupling (p.4j ) the lowest order term does not 
involve derivatives: 

£*V = ^(^ 2 7r7r + / 1 )(Tr#$p^ + Tr$#p^) (B.8) 
- Z U + f 2 ) Tr &pr R Q PLlt + Zf 3 (Tr $t$ - p ) Tr {p^p Lp + p R p Rp ) . 
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Our conventions are such that for /i = / 2 = /3 = the couplings ( |B.4|) , ( p.7| ), ( |B.8|) and an 
appropriate scalar kinetic term can be written in terms of a gauge covariant derivative 

= <9 M $ - -j=g p7r7r p RfJ ,<$> + -^g P nn$PLn (B.9) 

as ZTy (D^). For g P g q = g pwn this also extends to the couplings to quarks. 

Chiral symmetry breaking by the expectation value of $ leads to a mixing between p£ 
and p R . In the absence of quark masses the p-mass matrix reads 

~ { -2 fa, + h ) Zo% , ml + 2 + /,) Zol j ' (B ' 10) 
The mass eigenstates are the vector and pseudovector mesons (and associated scalars) 

Pv = -7z(Pr + Pl) 

Pa = 77= (Pb - Pl) ■ 

They transform under charge conjugation as 

c T 
Pv — > —Pv 

c 7 (R12) 

Pa — > P A 

and we conclude that the transversal parts of py and p^ A describe the 1 and 1 ++ octets 
and singlets of the light meson spectrum. The mass of the vector mesons is given by 

M v = m 2 p + 2{h- f 2 )Za 2 Q (B.13) 

whereas the pseudovector mass reads 

M\ = m 2 p + 2 (2g 2 p7Tn + /i + / 2 ) Za\ . (B. 14) 

We have put here a bar on m\ in order to indicate that the relation with physical axialvector 
masses involves an additional wave function renormalization 

M\ = ^. (B.15) 

In fact, chiral symmetry breaking induces also a difference in the kinetic terms for py and 
Pa by invariants of the type Tr &F r Q>Fl ij ,v with F^ V R = d^p u LR — d v p^ L R . Combining the 
most general term involving up to two powers of $ 

A£ ki n, P = ^ Tr &Fg*F Lia , + ^ Tr (&<S>F£"F L , U + M^iV) (B.16) 
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with ( p.2|) one finds for the kinetic term relevant for the spin-one bosons in the limit of 
vanishing quark masses 



= 4r Tr F v F v^ + ^ Tr Ff iV (B.17) 

o o 

with 

= Z p + {(32 + 13^1 

Z A = Zp+ifc-^al. (B.18) 

The wave function renormalization Z p can be fixed by convention. If we choose the nor- 
malization of pi r such that the kinetic term for p v has the standard form {Zy = 1) there 
remains nevertheless an additional parameter Z A multiplying the kinetic term of p^ which 
typically differs from one. Similar effects can influence the effective kinetic terms for the 
spin-zero components ~ d^p^. 

We observe that for g 2 + / 2 > the vector octet is indeed lighter than the pseudovector 
octet. For \f 2 \ <C g pn7T the mass splitting can be related to the pirn coupling and therefore 
to the p lifetime by 

M\ - M v = 4 (g%„ + f 2 ) Zal = ~ (2f K + J n f g% n ^-x p (B.19) 

with 

x p = 1 + . (B.20) 
The mass splitting within the octet because of nonvanishing quark masses can also be un- 



derstood from the interactions ( |B.8|) and ( |B.16| ) inserting $ = a + |(u>A 3 — y/3v\ 8 )Z h 1 ^ 2 . 



(This can be used to determine the parameters appearing in these expressions.) 

It is instructive to write the cubic coupling (|B.4|) in terms of mass eigenstates by using 



£<£ 2 p = t:9 p ^Z Tr \ [d p %% - + <^$ s $ s - $ S <9 M $ S ] p v 



2 

+ i 



2 

d(i*&p*&s ^pd p ^s *&sdfj,*&p dfJ&s^p ~\~ g i,Xsd p Xp ^ ^XsXp) 
2 , 1 (B.21) 

+ -j= (Xsd^p - d^ Xs % - Xpd^s + d^Xp^s) 



+ V2g pn7T Za Tr j + ~j=Xp^j d^pA 



We note a mixing of the longitudinal component d^p^ with the pseudoscalar mesons $ p , Xp 
for <7 > 0. This is possible since d p p^ A represents a 0~ + state. In contrast, d p p v transforms 
as + ~ and a mixing with ++ states $ s is possible only for the charged scalars. It is induced 
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by a nonvanishing expectation value ($ s ). As required by C and V invariance the vector 
mesons have cubic couplings to two pseudoscalars only if those are distinct, e.g. there is a 
p°7r + 7r~ but no p°7r 7r° coupling. Typical decays of pseudovectors involve the coupling of p A 
to one pseudoscalar and one scalar. Additional cubic couplings involving two (pseudo)vectors 
and one (pseudo) scalar are generated by ( |B.8| ) if a is inserted for one of the fields $. 

For an estimate of the effective ^-interactions induced by the exchange of p-fields we 
have to solve the field equations for p v and p A as functions of $. The result is then inserted 
into the action thus eliminating p M and replacing it by functions of $. We will do so keeping 
only terms linear and quadratic in $. This is sufficient for invariants containing up to four 
powers of $. For the contributions from the quartic term ( |B.8| ) we keep only the lowest 
order $ = diag(cfo), i.e. we neglect the mass splitting within the octets. The resulting field 
equations for p v and p A are 



d u dj; + ll - — J 8*8, + M v 5°) p a v b a = -ig^zfa&Q - &dfi 



ab 



-d»dj; + ( 1 - -J ^ + M 2 A s;j Pa \ = i 9pmt Z{d^ - 



(B.22) 



If we omit first effects from chiral symmetry breaking one has My = M A = m 2 . Inserting 
( |B.22|) into ( |B.2|) and ( |B.4|) and keeping only terms involving up to two derivatives, the 
vector and pseudovector mesons contribute to the kinetic term ( f4.6| ) only a structure 

An 2 7 2 

x (p) = _ y P n, (B 23) 

v m 2 

For <t > the SUl(3) X SUr(3) symmetry is spontaneously broken. The effective 
interactions mediated by py and p A still preserve the vector-like SU(3) symmetry if quark 
mass effects are neglected. In this approximation it is most convenient to give directly the 
contribution from p-exchange to the wave function renormalization constants Z m , Z p , Zh 
and Z s as well as u> m , u pm and u>h- They can be read off from (P-21| ) 

Aq 2 . 



Z ip) 

m 


= & 


7 {p) 


= W 




_ *g, 



— =2~~ Z a 

= 4 p) = o 

U Z " - (B.24) 
cro 



m a Z£Z„. 

tn) ^ SpTTTT Z 2 _ 

VQM 2 A (Z h Z m Z p )^ 
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2 

This yields the same expressions as inserting (|B.23|) into ( 8.10|) if m 2 is replaced by M A . 
Indeed, to linear order in the quark masses there is no contribution to the wave function 
renormalization constants from the exchange of the vector field p v . Only the exchange of 
the spin-zero component d^p'X induces the corrections ( |B.24j ). This can easily be understood 
from the structure of the interactions (p.21| ). Possible contributions to the quadratic terms 
defining Z m , Z h , oj m etc. can only arise through terms in the field equations for p which are 
linear in $ — ($). By Lorentz-invariance such terms must be ~ / cftxd^p 1 * ~ / d A x(§ — 
($))<9 M p M . The discrete symmetries C and V allow to this order only a term ~ Qpd^p^ which 
corresponds to the mixing of the _+ states and a structure ~ [$ s — , ($ s )] d^p v for 
the mixing in the scalar sector. The mixing in the pseudoscalar sector receives contributions 
~ o"o, v and therefore contributes to the terms in ( |B.24| ). In contrast, the mixing in the scalar 
sector vanishes for v — and gives therefore only a correction ( 11.2 ) to the scalar kinetic 
terms which is quadratic in v 



1 7 2 n 2 

The contribution to Z&) (|B.24|) is known as the "partial Higgs effect" . Inserting (|B.14|) 

2 

for M A one finds for f\ = fa = and Z = Z^ 

2 

Z m = Z v + Zg> = Z^ (B.26) 

M A 

and observes that Z m vanishes for m 2 — > 0. In this limit the symmetry SUl{3) X SUr(3) 
becomes an exact local gauge symmetry. As a result of the Higgs effect the pseudoscalar 
octet disappears from the spectrum. What remains are massive pseudovector mesons which 
acquire their mass through spontaneous symmetry breaking of the axial SUa{3) by Wq ^ 0. 
In the real world, however, the local gauge symmetry is explicitly broken by the mass term 
m 2 and by the deviation of fi and fi from zero. Below we will take Z different from Z^ 
and this again violates local SUl(3) x SUr(3) symmetry and modifies ( |B.26| ). The partial 
Higgs effect corresponds to the mixing between the two h octet states contained in $ p and 
dnP A . We may further improve the estimate (|B.24|) by taking into account the contributions 
~ q 2 on the left hand side of (p.22| ). The propagator of the divergence of p A differs from the 
one for the pseudovector mesons, since the kinetic term is given by the "gauge fixing" term 
~ ^- in ( |B.2|) . (There are similar differences in higher order in the momentum q 2 and in 

contributions to the kinetic terms from chiral symmetry breaking.) The expression for M A 
appropriate in ( |B.24j ) should involve the inverse propagator for the longitudinal component 
of p A , Zpq 2 + M A , where to lowest order Zp = a^ 1 . It should be taken at a momentum 
scale corresponding to the light pseudoscalar octet masses, q 2 = —m 2 ^- We also define a 
renormalized mass parameter 

M 2 P = M A /Z P . (B.27) 
Inserting this into ( |B.24|) yields the relation 

prfute&K+j*) = -\ ^/ p z 2 ( 2 ^+7,) 2 # 

Z mJ v ' 3 M A - Z P m 2 m K ' 
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6 Z M 2 A -M V /Z A 
x p Z m M\ - m 2 m Z P /Z A 



(B.28) 



where we used (|B.19|) for the last equation. We note that unless Zp/Za = M\/M P is very 
large the precise value of Zp has only little influence. In fact, we will see that Zp/Za ~ 0.35 
(cf. table The average masses of the vectors and pseudovectors are given by (M p = 
770 MeV, M K . = 892 MeV, M 0l = 1230 MeV, M Kl = 1340 MeV) 



My 



-ML + -Ml - 
3 K 3 p 



(853 MeV) 2 
(1300 MeV) 2 . 

This leads to a quantitative estimate for the p-exchange contribution to u m v 



M 2 A 



-Ml +-M 2 - 
3 Kl 3 ai 



(B.29) 



4 Z Uk ~ U (1 - M v /Ml) 
x p Z m (2f K + f n ) (1 - ml/ Ml) 



My 



M\ 



My 



Z A 1 



0.28 Z 



x p Z n 



0.75 ( — 

-Za 



(B.30) 



where in the last step we have neglected the weak dependence on Mp. At this stage we see 
already that a typical order of magnitude for u)$v is around —0.2 which is quite close to 
what is needed to explain f v and (see sect. ||). 

For a more detailed estimate we have to determine the appropriate choice of Z. After 
elimination of pi by solving its fields equation the relevant terms in the effective interaction 
for $ and py are 



C 



Z m Tr(d^d^ 



ZTt 



As a result of the partial Higgs effect we notice ( B.24|) the appearance of 



^9 pirn ^-r2 



07 



' P ™M A 



Zpm, 



(B.31) 



(B.32) 



instead of0 Z v in front of the kinetic term for $. Within the lowest order derivative approx- 
imation employed here we therefore have to identify 

Z = Z m (B.33) 

in order to get the standard coupling of py to the renormalized scalar field. Indeed, the 
lowest order pw interaction now takes the form 



(B.34) 



2:S 



More precisely, stands in this appendix for Z v + U v ao + X v a\ where X„ is the part of X v which 
is unrelated to "p-exchange" . 
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Xp 


z v 


Z A 


m a 
MeV 


Z P 


Mp 
MeV 


MO 

MeV 


( n \ 

v 






1.0 


1.00 


0.67 


1068 


0.00 






-0.14 


0.73 


0.00 


0.8 


0.99 


0.62 


1030 


0.22 


2196 


1852 


-0.17 


0.71 


0.08 


1.0 


0.99 


0.67 


1069 


0.22 


2278 


1944 


-0.15 


0.73 


0.07 


1.2 


0.99 


0.72 


1106 


0.22 


2358 


2031 


-0.14 


0.74 


0.07 


1.0 


0.9 


0.67 


1069 


0.16 


2673 


2282 


-0.15 


0.73 


0.05 



Table 9: The table gives the "leading mixing" results for various parameters related to 7 7-7/ m ixing. Only 
in the last line we use a fixed value Z p /Z m = 0.9 corresponding to a nonvanishing U v in (S.10). 



as can be seen by inserting p v = p^Ti, $ = l^ 1 / 2 /^ exp {j^f^), hj,k = 1,2,3. With the 
choice ( |B.33| ) the interactions ( |B.4| ), ( |B.7| ) can be combined for fx = f 2 = into a covariant 



kinetic term for $ 



4, y = Z m Tr(D^y{D^) 

Dp<S> = dp<5> - l -g p ^ [p V p, $] . (B.35) 

We observe in passing that because of the violation of gauge symmetry by the mass term 
~ m 2 p pAPAfi one cannot have simultaneously local gauge invariance of the interaction terms 
with respect to SUl(3) x SUr(3) (corresponding to Z = Z v ) and the vector-like subgroup 
SUv(3) after spontaneous symmetry breaking (corresponding to Z — Z m ). We could, of 
course, replace in (|B.4p g p7T7r Z m = g^Z^. The quartic interactions ( |B.8j ), however, would 



then be proportional to Z m g 2 p7T1T = Z^g 2 ^ (l H — v z m ). In this normalization (p.35| ) cor- 



responds to nonvanishing f x = f 2 = v z m and a different normalization of the gauge 
coupling. We will consider here the case of approximate SUy(3) gauge symmetry ( |B.35| ) in 
contrast to part of the literature where the choice Z = Z v is adopted while f\ and f 2 are 
neglected. In our case, the interaction terms are only invariant with respect to global chiral 
SUl(3) x SUr(3) transformations. Local S77y(3) invariance arises effectively only after chiral 
symmetry breaking if the symmetry breaking terms m 2 , fi, etc. are neglected. 

The wave function renormalization Z& can be evaluated from the average mass of the 



pseudovector mesons ( |B.19| ) 



M v + 4^ct 2 x 
M 2 



Z A = v u " . (B.36) 



If we neglect the weak dependence of tuffv on Z P , ( |B.30| ), ( |B.33| ) and ( |B.36| ) are already 



sufficient to compute uj^v as a function of x p . For a more accurate estimate we have to 
determine Zp or Mp which will be done in appendix |C] (cf. ( |C5| )). This allows to determine 
Z A , M A , Z Pl M P and 

LtJ m >v for given x p . The results are displayed in table |9] for different 
values of x p and M (cf. appendix |C|) . The first five lines are evaluated in the leading mixing 
approximation as described in the end of sect. 0. The last line assumes in addition U v 7^ 
in flS.lOD and corresponds to line (d) in sect. 13|. The main uncertainties in uj^'v arise from 
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Xp 7^ 1 and from possible higher derivative contributions to the decay p — > nix which would 
modify the value of g p7rn . In addition, there are also corrections from the mass splitting in the 
vector- and pseudovector octets. The latter are suppressed by an additional power of quark 
masses. We believe that the values for Z m jZ^ in table || give a rather realistic estimate for 
the quantitative importance of the "partial Higgs effect" . For the values in table [| we have 
assumed uj m = u>!$ for the computation of f K — f % , cr , etc. If we use instead the values of 
line (b) in table |, typically decreases by 0.02 thus reaching values close to —0.2. In 

conclusion, typical values for are in the range between —0.14 and —0.17. This coincides 
more or less with the values needed for an explanation of the observed decay constants f v 
and f n A We also note that the partial Higgs effect reflected in the deviation of Z m jZ^ from 
one leads typically to a 30% correction in the kinetic terms. 

Let us finally turn to the size of the mixing effects in the scalar sector parameterized by 
( h ( jggp . Despite their suppression by a factor ~ v they may be quantitatively relevant 
because the propagator of the exchanged quantum now involves the mass of the vector instead 
of the axialvector mesons. Since this mixing concerns only the Kq mesons (cf. ([11.2 )) we 



use for a quantitative estimate of (h the mass M K * = 892 MeV and replace in ( |B.25| ) M v by 
M|-« + Zpql = M\* — ZpMl*. (Here Zp differs from Zp by contribution ~ <7q similar to 
( |B.16|) ). This implies 

a P )2 _ i z m ( g P U7 K -L) \ 2 L ~ 7 M k \ 1 „ Q.Q185 

ChV -9X{ V'^MfJ ~ 1 - 2.57Z P (R37) 



and we note that this formula is not applicable for Zp around 0.4. There M\.jZp is of the 
same size as M\, and the mixing of states cannot be described by a derivative expansion 
anymore (cf. appendix 0). Also the imaginary part of the propagators of unstable particles 
have to be taken into account. Because of the uncertainty in the value of Zp it is difficult 
to give a precise quantitative estimate of C^v 2 . We only know that this quantity is positive 
and exceeds the value for Z P = 0. This leads to the estimate (h^ 2 ~ 0.02. 

We conclude that the "classical" exchange of spin one mesons only contributes to the 
effective quartic (and higher) interactions of the pseudoscalars but does not modify their 
effective propagators. They are therefore not relevant for the investigations of the present 
work which concentrate on masses and mixings. The only contributions from higher states 
to the propagators concern the mixing with higher 0~ + states for $ p and \ v and additional 
scalar states for $ s . Such states are contained in <9 M p^ and d p p v . These mixings contribute 
to quantities like X~, Zh/Z m etc. A brief general discussion of mixing effects is given in 
appendix y. 



C Mixing with other states 

In QCD the pseudoscalar and scalar mesons described by the field $ are only part of a 
rich spectrum of quark-antiquark states plus glueballs and possibly also qqqq states. There 
are strong couplings between the various states and their physics therefore influences the 
behavior of the h and ++ particles described in this work. On the level of effective 
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propagators which are the main subject of this paper we have to consider mixing effects 
with other CT + or ++ states in the spectrum. For the pseudoscalar octet this concerns 
only the higher mass _+ octets, whereas for the scalar octet we also have to consider the 
possible mixing with qqqq states in the ++ channel. We mention that it is not crucial in this 
context if the two-meson or four-quark states (with qqqq quantum numbers) correspond to 
"particles" like the a (980), / (980) or not — the composite fields describing the qqqq states 
may also have "propagators" without a pole. Finally, for the pseudoscalar rf we also have 
to include a mixing with pseudoscalar glueballs. We collectively denote these additional 
resonances as "higher states" . 

The general method for dealing with the higher states is to integrate them out and 
to compute an effective theory for $ alone. The investigations of this paper should be 
understood in this context. There are various methods for integrating out the higher states. 
One consists in computing first the effective action including additional fields for the higher 
states. In a second step the field equations for these states are solved for arbitrary values 
of $. The resulting "classical fields" are functionals of $ and can be reinserted into the 
effective action, thus leading to an effective action which depends only on $. The discussion 
in appendix [B] can serve as an example. For those results of the present paper which are only 
based on symmetries it is actually not necessary to perform the integration of additional fields 
in practice. Nevertheless, some insight in the origin of some of the constants of the effective 
action, like X~ V v etc., can be gained by considering the possible form of the effective action 
including additional states. We should point out that we neglect throughout the imaginary 
part of the two-point functions which is due to the decay of unstable resonances. This 
approximation may become invalid in the immediate vicinity of poles in the propagators. 

We have already encountered the mixing of h states in the discussion of the longitudinal 
component of p A in appendix |B[ Let us rephrase this with a somewhat different perspective 
by introducing a field 

T P = d v p v A - l ~Trd v p v A (C.l) 

for the additional h state. With this normalization the propagator for Tp can be approx- 
imated by q 2 Gp(q 2 ) with Gp 1 = Zpq 2 + M A (cf. appendix |B|). The inverse propagator for 
the coupled system of \/2& p and r P contains off-diagonal terms 

, G-\q) , G-p(q) 
G-\q) = V V „ (C.2) 

{ Gpl(q) , Gp 1 (q /g 2 1 



which are responsible for the mixing. From ( B.21 ) one finds 



G- P {q) = G P l(q) = 2g p7Tn Z m a . (C.3) 

A similar mixing occurs between \p an d the singlet t' p contained in Trd^p A . It is no accident 
that the mixing vanishes for a = or q 2 = 0: In the limit of unbroken chiral symmetry 
((To = 0) the fields <3> and p belong to different representations of SUl(3) xSUr(3) and cannot 
mix. Also for a ^ and vanishing quark masses $ p describes Goldstone bosons which can 
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a ,v 



00 



m rp p m t p p 

Figure 9: Feynman diagrams contributing to 77-77' mixing due to exchange of <9 M p^. 



only have derivative couplings. By construction the quark mass terms only appear as source 
terms for $. 

It is equivalent to diagonalize the matrix ( |C.2| ) or to eliminate rp by solving its field 
equations for rp[$] which is more adapted to our purpose. The elimination of rp gives an 
additional contribution to the effective inverse propagator G~ 1 (q) + AG~ 1 (g), namely 



AG-/(q) = -G- p (q)G P (q)q 2 



(C.4) 



There are also contributions to the off-diagonal m — p kinetic term related to 77-7/ mixing 
which are represented graphically in fig. |^. From an investigation of G^ 1 + AG" 1 we can 
obtain some general insight in the structure of mixing effects. First, we observe that for 



G^ 1 



Z^q 2 + M the effective propagator G = (G 1 



AG- 1 ) 



1 typically has two poles 



corresponding to the values of q 2 for which the determinant of G~ x (q) vanishes. In the 
vicinity of the lower mass pole G^ 1 + AG^ 1 can be approximated by a typical one particle 
two point function. The propagator vanishes at the value q 2 = —Mp where Gp 1 (q) has a 
zero. One observes that due to the particular factor of (q 2 )~ l in the inverse rp propagator 
( p.2|) the value of Mp is always larger than both masses corresponding to the location of 
the polesQ. In addition, the residue of G at the pole with the higher value of — q 2 has the 
opposite sign as for the lower mass pole. (The higher mass pole does not correspond to a 
stable particle even within our approximations.) If we denote by q 2 = —Mq the location of 



the higher pole one finds with (|C3[) and M 



M 2 P 



M\/Z t 



Ml 



1 + ^i Z t[ M Z 



Z,, 



(C.5) 



An estimate of Mq is not obvious and subject to large uncertainties. One can then determine 
Z^/Zm from (|B.32|) , (B.33| ) and solve the resulting system of equations in dependence on 



Zp. Results are displayed in table || for values of Zp in a range for which M n comes out 
with a reasonable size in the leading mixing approximation (fig. |8]). We note that much 
larger values of Zp lead to a very large mixing in the 77—77' sector and completely destroy any 
reasonable picture. We also show in table three different values of x, 



24 Without the (q 2 ) 1 factor Mp would be inbetween the two poles. 
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Second, it is clear that for real octet or singlet fields the off-diagonal elements must 
be real and equal. Integrating out the additional fields gives a negative contribution to the 
coefficient ~ q 2 in the quadratic term for <&—($) as long as Gp(q) stays positive. The mixing 
gives therefore a negative contribution to Z m (|8.10|) and explains why X~ is negative. Third, 
from ( |C.4j ) we learn that the mixing effects are proportional to the propagator Gp. This 
suggests that mixing effects with light additional states are particularly important. Fourth, 
the mixing effects also contribute to higher derivative terms in the effective action for $. 
Expanding Gp(q) around q 2 . gives 



AG 



-i 



-4 



„2 V2 ^2 

•3 pTTTV Tfl 2 

M\ + Z P ql 



Zp {q 2 ~ qg) 
M\ + Zpq 2 



(C.6) 



where the first term corresponds to Z^q 2 as given by 
using (p.24|) we find a contribution 



Comparing with ( p.l2| ) and 



-Trip) 



4 m f K ~L M 2 P -ml 



(C.7) 



The momentum dependence from the effective propagator Gp may be particularly im- 
portant for the 7/ since its mass is closest to Mp. The resulting g 2 -dependence of the 
quantities appearing in ( |6.26D has therefore to be treated with care. We remark that typical 
masses in the neutral sector are higher than in the flavored one and therefore guess Mp 
around 2000 MeV with a large uncertainty. In terms of the unrenormalized fields the inverse 
propagator in the flavor neutral sector takes on the form 



Z^ + X-{q 2 )al] q 2 + M 2 Z p 



u(q 2 )q 2 



zll 2 z l J 2 



Z v + X-(q 2 )a 2 } q 2 + hll 



(C.8) 



Here we assume that the nontrivial momentum dependence beyond the approximation linear 
in q 2 arises dominantly from propagator effects contained in X~(q 2 ) and Cj(q 2 ) according to 



X-{q 2 ) 
u(q 2 ) 

Our conventions imply 



{M 2 



(M 2 + q 2 ) 
M 2 )Uq 2 ) , UQ 2 ) 



M 2 



M 2 



Mp + q 2 



(C.9) 



-Mi 



x: 



and we remind that we have defined both, 
approximation implies 



Zp and Z 8 , 



for q 2 



-M 2 . 



(CIO) 

The leading mixing 



x: 



Z m Z s 



x: 





M 2 



—2 

r 1 J a 



(G.U) 
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With this definition the only higher derivative effect in the matrix ( |6.26| ) for q 2 = —M 2 
appears in the factor f~ l (-m 2 m ) = (Mp - m 2 m )/(M 2 P - M 2 ) = 1.02, ( $M ). Here we have 
used the estimate Mp 
unity is small 



2670 MeV from table |^. Also the resulting deviation of Z p /Z m from 



Z„ + x- 



-M 2 )a 2 



Zm + X~ 



1 + 



x: 



1 (2f K + f n )(M 2 -m 2 J 



fn) ( M P ~ M v) 



-M 2 ) - X- 



1 + 0.05uo m v 



(C.12) 



and compatible with the linearization for Z s according to ( C.7|) . (The numerical value for 
= —0.20 is Z p /Z m = 0.99.) There is no modification of the decay constants f v& and /^ , 
( |A.24| ), (|A.25 ). On the other hand, the propagator effect in X~(q 2 ) and Cj(q 2 ) could lead to 



substantial effects for q 2 = —M 2 , depending on the value of M P : In the diagonal elements 



of the inverse propagator ( |6.26|) one has to insert 

Z^X-{q 2 )a 2 



z v + x-> 



-M 2 )a 2 



_ I, , .. (2/ (Q 2 +m 2 ) 



(C.13) 



whereas u is replaced by u)(q 2 ) ( |6.33| ), flC.9| ). Furthermore, the correct definition of the decay 
constants f^o and f v >s involves now the factor ( |A.26| ) 



(C.14) 



For the scalar octet an interesting possibility of mixing concerns states in the two-meson 
channels. In fact, the four-point function for $ p may develop resonance-like structures in 
the momentum range corresponding to the sum of two pseudoscalar meson masses. (These 
momentum dependent structures are not accounted for by the four-point function at zero 
external momenta described by the effective potential U .) Such resonance structures can 
be replaced by effective interactions with a composite ++ field ts- The effective two-point 
functions for t$ obtained in this way do not necessarily correspond to a propagating particle 
or resonance, since their real part may be strictly positive and bounded for all values of q 2 
on the real axis. Consider a generic structure for the mixing between v^2$s and t$ 



G-\q) 



Zq 2 + M 
b(q 2 ) 



b(q 2 ) 
c{q 2 ) 



(C.15) 



If c(q 2 ) has a zero for \/—q' 2 in the vicinity of the sum of two pseudoscalar masses one finds 
two values of q 2 for which an eigenvalue of G~ l vanishes. In the case of the isospin triplet 
they could be associated with a (980) and a (1320). On the other hand, the two-particle 
threshold could also be reflected by a finite enhancement of c(g) -1 or b(q) without a zero of 
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c(q). If c(q 2 ) dips in this momentum region to values smaller than b 2 (q 2 )/(M + Zq 2 ) the 
zero eigenvalue of G -1 will occur precisely in the threshold region, namely for q 2 = — Mq 
as determined by c(-Mq) = 6 2 (-Mq)(M - ZM^y 1 . After solving for r 5 [$ s ] the location 
of the single pole of [G" 1 ^) — 6 2 (g 2 )c _1 (g 2 )] -1 would then necessarily be found at — Mq in 

the threshold region. For M 2 jZ not too far from the two-particle threshold this effect could 
explain naturally why the isotriplet in $ s is found precisely at the 2K threshold! Mixing 
effects from AG^ 1 = —b 2 {q 2 )c^ 1 {q 2 ) are large in this case. Since this mechanism requires a 
critical strength for 6 2 c _1 not all members of the scalar octet have to be in the vicinity of 
two-particle thresholds. More precisely, the phenomenon of "threshold mass shifting" which 
induces mesons masses near a two-particle threshold occurs whenever M jZ is above the 
threshold and the quantity b 2 (q)/(c(q)Z(q)) makes a strong enough jump in the threshold 
region. An alternative way of looking at this "threshold mass shifting" notes that the loop 
contribution to the two-point function for the do becomes important if the mass is close to 
the sum of the masses of the two pseudoscalars circulating in the loop. 

Even without a detailed discussion of the complicated analyticity properties we conclude 
that for both alternatives the effective inverse propagator for the isotriplet in $ s should have 
a zero at the observed ao(980) resonance. In this momentum region the mixing effects with 
two-kaon states are expected to be very strong. No detailed understanding of the properties 
of the ao(980) seems possible without incorporating the two-kaon channel. In addition, the 
effective inverse propagator may (or may not) have a second zero corresponding to the possi- 
ble resonance ao(1320). In this momentum region the mixing effects should be much smaller 
because of the larger value of Zq 2 + M — c(q). We note that both resonances are described 

2 

by the same value o^| M , but different effective Z h and oo h . The two different associations 
a (980) vs. a (1320) in the main text become in this case only two facets of the same story. 
If the a (1320) exists the values characterizing the potential, like m\, A 2 , A 3 , etc. should 
be independent of the identification of the isotriplet. The actual differences in these values 
are then a measure for the influence of neglected terms. We find that these differences can 
indeed be small if the mixing is large enough for the ao(980). This is compatible with the 
existence of the ao(1320) real resonance. 

Finally, we turn to the mixing of a pseudoscalar h glueball g with the rf. This is 
particularly interesting in view of a possible experimental detection of g. Because of the 
anomaly the rf is not a Goldstone boson for vanishing quark masses and the off-diagonal 
element in the mixing matrix may not vanish for zero momentum. We introduce for the 
glueball a pseudoscalar singlet field g with an effective action 



The coupling between $ and g conserves all symmetries (V(g) = —g, C(g) = g). The real 
coupling h g i may depend on the momentum of g. Expanding h g i(q) around q 2 = —rr? m the 

25 It is convenient to choose composite fields such that c(0) = 0. 



£g = ^gd^g + -m 2 gl g 2 + h gl ug 



10 




(C.16) 
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mixing with the glueball contributes to the effective potential for $ (cf. 



2 91 



CO 



whereas for the kinetic term (03) it induces a coupling 



Vi 9) = h 2 (q ) \m 2 gl + q 2 



dfl^i r 2 9 

g^-Wo) [m gl + q 



(C.17) 



(C.18) 



It is probably difficult to disentangle (p.!7f) from other contributions to the potential. On 
the other hand, a determination of the size of the parameter V v will put restrictions on h g i. 
Since h g i is directly related to the mixing between rj' and g one may obtain from it interesting 
information on the decay of the pseudoscalar glueball into mesons or photons. To lowest 
order we can use , 

to = -V6a 2 ZPp (C.19) 

and obtain the inverse propagator for the r/-glueball system as 



G-\q) 



q 2 + ml 

2 a O a gl 



q 2 + m 2 gl 



(C.20) 



For not too large h g i the mixing angle between g and rj' is suppressed by the small ratio 
Vo/m g i 

2 

9 l /° 2 ^10' 3 ^- (C21) 
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